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Preface 



It has been felt for a while that our graduate students do not get an opportunity to get an 
exposure to the non-perturbative, background independent quantum theory of gravity at a 
pedagogic level. Although there are several excellent reviews and lecture notes available, 
an opportunity for complementing lectures by discussions is always an added bonus for the 
students. With this in mind and taking into account of the background preparation of the 
students, the School on Loop Quantum Gravity was organized at IMSc, for a period of 10 days. 
The first 5 days were devoted to the basics of connection formulation and loop quantization 
up to sketching steps involved in the quantization of the Hamiltonian constraint. The next 5 
days were devoted to the applications to quantum cosmology and to the black hole entropy. 
In all 20 lectures and 10 tutorials were planned, however some tutorials 'became' additional 
lectures. 

These notes are an expanded version of the topics that I covered. In particular the material 
of chapter 2, sections 4.2.1, 4.2.2, 4.2.3 and appendices 5.1, 5.3, 5.4 have been added. Email 
discussions on the sections of chapter 4 and appendix 5.4 with Abhay Ashtekar, Martin 
Bojowald and Madhavan Varadarajan have been very helpful and are gratefully acknowl- 
edged. There could still be some differences in the perceptions and formulations, what I 
have presented is my understanding of the issues. 

The other main lecturers at the school were: Prof Amit Ghosh, Saha Institute of Nuclear 
Physics, Kolkata; Dr Alok Laddha, Raman Research Institute, Bengaluru; Parthasarathi 
Majumdar, SINP, Kolkata. In addition, Prof Romesh Kaul, IMSc, Dr Kinjal Banerjee, 
IUCAA, Pune and Ayan Chatterjee, SINP, Kolkata also gave a few lectures. Amit and 
Alok discussed the connection formulation and loop quantization up to the basic steps in 
the quantization of the Hamiltonian constraint. Partha, Amit and Ayan discussed classical 
formulation of isolated horizons and entropy associated with them. Romesh discussed the 
possibility of a 'vacuum structure for gravity' and Alok also briefly discussed the Brown- 
Kuchar dust model. It is a pleasure to acknowledge their contributions. At least some notes 
of the topics covered by these will be available in not-too-distant a future. 

The funding for the school was provided by the Institute of Mathematical Sciences under 
the XI th Plan Project entitled Numerical Quantum Gravity and Cosmology. It is envisaged 
that a more specialized workshop at an advanced level will be held in the summer of 2010 
with the possibility of a similar School being repeated one more time. 
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Chapter 



1 



General 



Remarks: 



1.1 Why a 



Quantum Theory of Gravity? 



This is not a rhetorical question but it is intended to identify physical context in which the 
classical theory of gravity, specifically the Einstein Theory of Gravity also referred to as 
General Relativity, is inadequate and calls for an extension. One has met with inadequacies 
of classical theories many times and has seen how their quantum versions have alleviated 
the inadequacies. For example, the classical theory of charges and electromagnetic fields was 
quite adequate until the hydrogen atom was found to have a central nucleus with an electron 
going around it. Classical theory predicts that since the electron is necessarily accelerated, 
it must radiate away its energy and spiral into the proton. Indeed in about 10~ 9 seconds (!) 
the classical trajectory of a (bound) electron would 'end' in the proton. We all know that 
this is physically wrong and the atoms are known to be stable for billions of years. We also 
know that the 'fault' lies not with the 'Coulomb law' (which does get modified) but with 
the classical framework of using well defined trajectories to describe the dynamical evolution 
for both the electron and the electromagnetic field. Inadequacies of classical theories are 
also revealed in the black body spectrum, specific heat of solids at low temperatures etc etc 
and the appropriate quantum theory of matter and electromagnet ism cures these problems 
i.e. gives results consistent with experimental observations. The quantum nature of other 
interactions such as the strong and the weak is also verified in nuclear and particle physics. 
What about the gravitational interactions? 

Gravitationally bound 'atoms' can also be considered. If gravity is described in the Newto- 
nian manner, there is no gravitational radiation from an accelerated motion and the inward 
spiralling problem will not arise. But Einstein's theory of gravity is very different and ac- 
celerated sources do radiate away energy and the stability issue re-surfaces. Of course the 
'weakness' of gravitational interaction does not threaten the existence of such gravitation- 
ally bound atoms if the decay time is larger than the age of the universe, but in principle 
possibility exists. In fact Einstein did suggest a need for a quantum theory of gravity [1] 
almost immediately after GR was constructed^- 

General relativity however uncovered two distinct contexts in which the theory calls for an 

1 At that time, the universe was supposed to be eternal and hence, however small the gravitational 
instability, existence of stable atoms would threaten GR unless gravitational radiation is also terminated at 
certain stage. 
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extension - the context of (i) cosmological and black hole singularities and (ii) entropy of 
black hole horizons. Let us take a little closer look at these contexts. 

The cosmological context: Under the assumption of homogeneity and isotropy, the space- 
time metric is described in terms of a single dynamical variable - the scale factor. As long 
as the energy density is positive and the pressure is not too negative (which is true for the 
properties of known matter), in an expanding universe (an observational fact), the scale factor 
vanishes at a finite time in the past. The universe has 'beginning', a finite age and the space- 
time curvature (or gravity) is infinitely large. Thus, a homogeneous, isotropic universe has 
singular beginning. If one relaxes isotropy but retains homogeneity, one has several types of 
space-times. There are now a maximum of six dynamical variables. The simplest of these, the 
vacuum Bianchi I space-time, has three 'scale factors' whose time dependence is given by the 
(exact) Kasner solution. This is also singular. As one evolves back in time, two of the scale 
factors vanish while the third one diverges. The most complex of these models, the vacuum 
Bianchi IX space, is also singular and the backward evolution has an oscillatory behaviour. 
Like the Kasner solution, two scale factors begin decreasing and the third one increasing. 
But after a while, the three scale factors change their behaviour and a different pair begins 
decreasing. This continues ad infinitum. If the non-diagonal metric components are included, 
then the directions along which contraction/expansion takes place are also 'rotated'. If one 
relaxes homogeneity as well, then a beautiful analysis done by Belinskii-Khalatnikov-Lifschitz 
(BKL), shows that there exists a general solution of the vacuum Einstein equations which can 
be described as smaller and smaller portions of the spatial slice behaving as a homogeneous, 
Bianchi IX solution. The BKL analysis in particular shows that singular solutions found in 
the simpler situations are not due to high degree of symmetry (homogeneity and isotropy), 
but even without such symmetries, there exist general solutions which are singular (diverging 
curvatures) and the nature of singularities can be extremely complicated. 

During the sixties Geroch-Penrose-Hawking used another approach to establish the Singu- 
larity theorems identifying conditions under which singularities are inevitable consequence 
of classical GR. For these theorems, singular space-times were defined as those in- extendible 
space-times which admit at least one causal (time-like or null) geodesic which is incomplete. 
Here incompleteness means that the geodesic cannot be defined for all real values of an affine 
parameter. There are three types of inputs in these theorems: (a) One restricts to a class 
of space-times which are causally well-behaved eg are free from closed causal curves. The 
space-times which are free of all causal pathologies and are fully deterministic are the so- 
called globally hyperbolic space-times, (b) The space-times are solutions of Einstein equation 
with the matter stress tensor satisfying suitable energy condition(s). This incorporates that 
idea that gravity is attractive (for positive mass/energy). These two types of conditions are 
general requirements for a space-time model to be physically relevant, (c) the third input is 
a condition that distinguishes specific physical context such as an everywhere expanding uni- 
verse or a gravitational collapse which has proceeded far enough to develop trapped surfaces. 
The presence of the last condition(s) shows that not every solution satisfying the first two 
conditions is singular (e.g. the Minkowski space-time). Thus, singularity theorems do not 
imply that gravitational interactions always produce singularities - the (c) type condition 
is necessary. While inclusion of (c) will imply singular space-times, it is not automatic that 
this condition is realized in the physical world. In our physical world however universe is 
expanding and it is widely believed that black holes also exist and hence condition (c) is 
realized in nature. Thus, physical contexts exist wherein classical GR is inadequate. 

Remark 1: The global hyperbolicity condition implies that the space-time is stably causal i.e. 
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a global time function exists such that each hyper-surface of constant 'time' is space-like. 
Furthermore, a time function can be chosen such that the space-like hypersurfaces are Cauchy 
surfaces. The topology of such space-times is necessarily K x E. A Hamiltonian formulation 
makes sense only in such space-times. Thus, not every solution of Einstein equation yields 
a physically acceptable space-time (i.e. causally well behaved or globally hyperbolic). A 
Hamiltonian evolution however constructs such space-times [E 

Remark 2: In specific solutions, one encounters singularities (regions of diverging curvatures) 
which are space-like (positive mass Schwarzschild solution, homogeneous cosmologies), time- 
like (negative mass Schwarzschild solution) or even null (some of the Weyl class of solutions). 
Singularities that arise in an evolution from non-singular initial conditions are the ones which 
strongly display inadequacy of the theory. Typically, these are the space-like singularities. 
Since the Hamiltonian formulation is an initial value formulation, it can "see" only such 
singularities. 

The Black hole context: Black holes are objects whose 'interiors' are inaccessible to far 
away observers. More precisely, these are space-time geometries that have a horizon which 
leave some regions out of bounds for asymptotic observers. The special class of stationary 
black holes are characterized by a few parameters - mass (M), angular momentum (J) and 
electric charge (say) (Q). Associated with their horizons are some characteristic parameters 
- area of the horizon (A), surface gravity at the horizon (k), angular velocity at the horizon 
(Q) and electromagnetic potential at the horizon ($). In the seventies, a remarkable set of 
"laws" governing processes involving black holes were discovered. If in a process a black 
hole is disturbed (by accreting mass, say) and returns to a stationary state again, then the 
changes in the parameters obey: 

SM = —5A + n H 5J + $ H 5Q , 8A>0. (1.1) 
8n 

These are very temptingly analogous to the laws of thermodynamics! Especially after one also 
proves that k is constant over the surface of the horizon. Bekenstein in fact suggested that 
area of the horizon be identified with the entropy of a thermodynamic system. This suggests 
that the surface gravity be identified with a temperature. Hawking subsequently showed 
that when possibilities of quantum instabilities are taken into account, a black hole can be 
thought of a black body with temperature T = and hence S = \ '■= GK. For all 
other systems we know that thermodynamics is a manifestation of an underlying statistical 
mechanics of a large number of microscopic degrees of freedom. What are these micro- 
constituents of the black holes? Notice that from far away, only the exterior of a horizon is 
accessible and so also parameters such as mass and angular momentum. All detailed memory 
of what collapsed to form the black hole is lost. So these micro-constituents must be distinct 
from the matter degrees of freedom. They must represent "atoms" of geometry. But classical 
geometry is continuous so how does a particulate nature arise? Perhaps, not just the specific 
dynamics given by Einstein equation is inadequate but the framework of classical geometry 
itself is inadequate. Note that black hole horizons are not regions of high curvatures and 
geodesic incompleteness occurs in their interiors. Thus, black hole thermodynamics is a 
qualitatively different situation. 

In summary, classical GR contains within its domain, physically realizable physical context 
where the theory is inadequate. At least one of its context involves highly dynamical ge- 

2 One could analytically extend such maximally Cauchy evolved solutions further. These will have Cauchy 
Horizons, an example being Reissner-Nordstrom space-time. 
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ometries with high curvatures, matter densities etc. Because of these features, it is hard to 
imagine how any perturbative approach can be developed in these contexts. Since gravity (or 
space-time geometry) is dynamical and a perturbative approach is unlikely to be suitable, it 
is necessary to have a quantum theory of gravity which of does not use any fixed background 
space-time in its construction. 



1.2 An Essential Feature of Classical Gravity 

Let us recall briefly that special relativity combines Newtonian notions of space and time 
into a single entity, the space-time (Minkowski space-time). The analysis of the geometry 
inferred in a rotating frame indicates that the geometry is non-Euclidean. Principle of 
equivalence, which identifies uniform acceleration with uniform gravity (in the Newtonian 
sense), then implies that gravity affects the space-time geometry and since matter affects 
gravity, it also affects geometry. In the final formulation of the relativistic theory of gravity, 
the space-time geometry, described by a metric tensor, is a dynamical (changeable) entity 
with the Newtonian gravity being a manifestation of the curvature. The law governing the 
matter-geometry interaction is encoded in the Einstein equation. That all observers are on 
equal footing to formulate the laws of physics implies that all quantities (and equations) be 
tensor fields (and equations) with respect to general coordinate transformations. Note that 
a general coordinate transformation corresponds to a change of chart in the framework of 
different iable manifolds. 

Such coordinate transformations however have another interpretation in terms of mapping 
of the manifold (or local regions thereof) into itself - the active diffeomorphisms. Under 
the action of such mappings, the pull-backs and push-forwards, generate "new" tensor fields 
from the old ones. That is, in a given neighbourhood, we will have the original tensor field 
and the one obtained via pull-back/push-forward. If x — > y(x) represents the mapping in 
terms of local coordinates, then the pulled-back (pushed-forward) quantities are related to 
the original ones in precisely the same manner as general coordinate transformation B if 
the dynamical equations are covariant with respect to general coordinate transformations 
(coordinate transform of a solution is also a solution), they must also be covariant with 
respect to the active diffeomorphisms i.e. a configuration and its transform under active 
diffeomorphism are both solutions if any one of them is. This has far reaching implications. 

The Einstein Hole Argument: The active diffeos can be chosen such that they map non- 
trivially only in some region ('sub-manifold') of the manifold. Choose a region which is 
free of any matter. Assume that the equations determining gravitational field and matter 
distribution are also tensor equations (i.e. generally covariant). Consider a solution which 
has certain curvature distribution inside our chosen 'hole'. Make a diffeo which is non-trivial 
only inside the hole and change the curvature distribution. This will also be a solution by 
covariance. Thus we get a situation that even though matter distribution is unchanged, in 
a region where there is no matter, we can have two 'different' gravitational fields i.e. matter 

3 Let 4> : p — > <j)(p) be a smooth map of a manifold onto itself. Given a function /' : M — > R, define 
another function / : M R as: f(p) := f'(4>(p)). This function is the pull-back of /' and is also denoted as 
/ := 4>*(f). Likewise, given a vector field X on M define a new vector field X' as: X'^f')]^) := X(f)\ p . 
The new vector field is called the push-forward of X and is also denoted as: X' — 4>*{X). Now introduce 
local coordinates x l around p and y l around <j>(p). It is easy to see that the components of tensors relative 
to these coordinates are related in exactly the same manner as though x — > y(x) is a change of chart. 
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distribution does not determine the gravitational field. But in the non-relativistic limit 
Newtonian gravity is determined by matter. So either the equations should not be covariant 
or in the absence of any matter available for 'marking' points of a manifold, the 'different' 
distributions of curvature must be regarded as describing the 'same ' gravitational field. It 
is the latter possibility that remains once the covariance of the equations is accepted. This 
in turn implies that it is the equivalence classes of solutions, with respect to the space-time 
diffeomorphisms, that correspond to physical reality. 

Note: The size of the hole in the hole argument is unimportant. Also, the metric description 
itself does not play a role; one could repeat the argument for any other field. All that is used 
is that the fields are tensors under general coordinate transformations (chart change), field 
equations are covariant and the fields are inhomogeneous within a hole. Values of individual 
fields at any, manifold point are irrelevant but values of fields at points specified by other 
fields are invariant and thus physical. 

Since 'points', not marked by any dynamical entity, have no physical meaning, the only, 
physically meaningful, questions are of relational nature. That is, it is physically meaningless 
to ask what is the curvature (or say electric field) "here and now" . The meaningful questions 
is what is the curvature where a certain field has a certain value. If we had any particular 
field to be fixed (non-dynamical), then with reference to that field we could ask the 'here and 
now' question. Such a field, constitutes a background. Note that the usual non-gravitational 
theories or in the perturbative treatment of gravity, the space-time geometry (metric) plays 
the role of a background. Since in the general relativistic theory all fields including the 
metric are fundamentally dynamical, such a theory is necessarily background free. The 
twin features of the framework namely all fields on a manifold being dynamical and the 
fundamental equation being generally covariant and deterministic, implies covariance with 
respect to active diffeomorphisms and physical characterization being in terms of 4-diffeo 
equivalence classes of fields Q 

This is an essential feature of general relativity, much more fundamental than the particular 
Einstein equations themselves. The challenge is to construct a quantum theory which faith- 
fully incorporates this feature i.e. a quantum theory of gravity must be background free (or 
at least recover background independence in the classical limit). 

This also poses challenges, because we have to construct observables which are space-time dif- 
feomorphism invariant. These alone could characterise specific equivalence classes of space- 
times and this problem is not understood even classically for spatially compact case and in 
absence of matter! Note that curvature invariants, although scalars, are local and not diffeo 
invariants. Hence these cannot be physical observables. Consequently identifying a physical 
state corresponding to (say) Minkowski space-time is much harder. 

For a more detailed discussion of these conceptual issues, see [2]. 



4 Since we have taken (differential) equations as specifying a presentation of a theory, the manifold cannot 
be thought of as a background, but rather part of the specification of the theory. If we take some transition 
amplitudes (among topological spaces, sets, ...) as specifying a theory, then the choice of a particular 
differential structure will constitute a background since it is also a dynamical entity. We will restrict to 
manifold category. This also explains why metric by itself is not essential for background independence of 
gravity, a dynamical tetrad with compatible spin connection would do just as well. 
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1.3 Towards the construction of a quantum theory 



Ultimately constructing a quantum theory of some phenomenon means specifying a state 
space - (projective) Hilbert space, identifying (self-adjoint) operators on it to correspond to 
physically observable quantities, a notion of evolution or dynamics such that in a suitable 
semi-classical approximation, the evolution of expectation values of a class of observables 
tracks the corresponding classical evolution with the quantum uncertainties less than the 
observational precision. Here the classical evolution is the one specified by the classical 
description of the phenomenon. Using such a framework, one can compute matrix elements 
of suitable observables or transition amplitudes etc. One familiar procedure is that of the 
canonical quantization. 

In canonical quantization, typically we have a classical phase space which is cotangent bundle, 
T*Q of some configuration manifold Q and the Hilbert space is the space complex valued 
functions of Q which are square integrable with respect to some suitable measure, d\x. When 
Q ~ M. N , we have the familiar L 2 (M. N , (i/^Lebesgue) which is unique thanks to the Stone- von 
Neumann theorem. This comfortable situation changes once Q becomes topologically non- 
trivial and/or becomes infinite dimensional. The former can arise due to constraints while 
the latter arises in a field theory. In relativistic field theory, the classical configuration space, 
Q, is (say) the space of suitably smooth tensor fields which is inadequate to describe the 
corresponding quantum fields which can be arbitrarily non-smooth. Typically, Q is extended 
to a quantum configuration space, Q, which should admit a suitable measure. 

For a quantum theory of gravity, there are two additional features - (i) we have a theory 
with first class constraints (i.e. a gauge theory) and (ii) we would like to have background 
independence. 

In presence of constraints, the quantization procedure is a two step process. In the first step 
one constructs a kinematical Hilbert space on which are defined the constraint operators. The 
second step aims to 'solve' the constraints to get a quantum theory corresponding to physical 
degrees of freedom. Again, typically, there are no vectors in the kinematical Hilbert space 
which are annihilated by the constraint operators and one is forced to consider distributional 
solution^. The space of distributional solution however is not a Hilbert space and another 
inner product needs to be defined on this space to make it into the physical Hilbert space. 
The choice of this inner product is limited by demanding that a suitable class of Dirac 
observables - operators which leave the space of solutions invariant - be self-adjoint. 

There are many choices to be made along the way. The requirement of background indepen- 
dence means that no non-dynamical fields should be used in any step. This poses a severe 
challenge to constructing even the kinematical Hilbert space. The connection formulation of 
gravity is of great help as the quantum configuration space of a gauge theory, A/Q - space 
of generalized connections modulo generalized gauge transformations - admits several mea- 
sures and the demand that the conjugate variables be represented by derivative operators 
essentially singles out a unique measure - the Ashtekar-Lewandowski measure, essentially 
constructed from the Haar measure on compact groups. One has a natural choice of Q 
and (non-unique) definitions of constraint operators so that the kinematical set-up is well 

5 Let 51 C i?kin be a dense subspace of the kinematical Hilbert space. Let 51* denote its algebraic dual 
(space of linear functions on SI) so that 51 C -ffkin C 51*. 51 is chosen so that it contains the domains of the 
constraint operators as well as of other operators of interest. Distributional solutions of constraints are those 
elements of 51* which evaluate to zero on all elements of 51 of the form C\ip),V if) E 51. 
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founded. We will see the details of these steps. 



We will begin with the Hamiltonian formulation of GR in terms of the metric (ADM for- 
mulation). Discover the redundant variables and make a canonical transformation to a set 
of new variables (connection formulation) which are amenable to background independent 
presentation. These will lead us to the holonomy-fTux variables and their Poisson bracket 
algebra whose representation theory will give us a Hilbert space. This will complete the 
first step in the construction of the quantum theory. Some of its novel features will be re- 
vealed through the properties of the geometrical operators. Our study of the basic formalism 
will conclude with the presentation of the constraint operators on the kinematical Hilbert 
space. The dynamical aspects will be studied through the simpler cases of homogeneous and 
isotropic cosmology leading to the Big Bang singularity resolution. The other application 
of quantum geometry, namely revealing the 'atoms' of geometry responsible for the black 
entropy will be discussed in the second week along with the loop quantum cosmology. 
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Chapter 2 



Classical Hamiltonian Formulation 



We are familiar from usual spacial relativistic field theories (say a scalar field) that solutions 
of the field equations can be viewed as an evolution of fields, their spatial derivatives and 
their velocities from one spatial slice to another one ("a Cauchy evolution"). In the general 
relativistic case, one has to deal with space-times other than the Minkowski space-time and 
eventually make the space-time itself to be 'dynamical'. However not all space-times support 
this notion of evolution. 

To have a well defined, causal (no propagation faster than speed of light in vacuum) and 
deterministic (given certain data at one instance, the future data is uniquely determined), 
the space-time must be free of causal pathologies such as (i) no closed causal (i.e. time- 
like or null) curves (excludes by chronology condition); (ii) no closed causal curves but 
causal curves which return arbitrarily close to themselves (excluded by strong causality); 
(iii) strong causality holds but when the space-time metric is made slightly 'wider', it is 
violated (excluded by stable causality). All such pathologies are absent in space-time which 
are stably causal i.e. admit a differentiable function such that d^f is a time-like vector field. 
This alone is still not sufficient to guarantee the possibility of a Cauchy Problem. For this, 
one needs Globally Hyperbolic Space-times. These are space-times which are which admit 
a spatial hyper-surface such that events to the future (past) are completely determined by 
data specified on it. Such space-time admit a globally defined 'time function' whose equal 
time surfaces are Cauchy surfaces. It follows further that such space-times must have the 
topology: M = R x E 3 . 

Thus, in order to have a well-posed causal theory of matter fields, the space-times must be 
globally hyperbolic. 

It is a non-trivial result of analysis of Einstein equation that Einstein equation can be cast 
in a Hamiltonian form such that if initial conditions are chosen to satisfy certain constraints, 
then corresponding Hamiltonian evolution generates a solution (space-time) of the Einstein 
equation^. 



lr This is local existence and uniqueness theorem for the Einstein equation. Since these are short time 
evolutions, one cannot guarantee that largest possible space-time constructed will be globally hyperbolic. 
However, if a globally hyperbolic solution is to exist, one can perform a time + space decomposition to put 
the equations in a Hamiltonian form. That Einstein equation admit a well-posed initial value problem is 
a necessary condition for globally hyperbolic solutions. That the equations are of Hamiltonian form is an 
additional, non-trivial property. This follows most directly via the Einstein-Hilbert action formulation. 
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Although a Hamiltonian formulations can be specified ab initio by giving a phase space which 
is a symplectic manifold, a Hamiltonian function and specifying evolution by the Hamilton's 
equations of motion, it is more common that a theory is specified in terms an action which 
is a functional defined on a set of fields on a space-time manifold. Typically this is expressed 
as an integral of a Lagrangian density made up of finite order derivatives of a set of tensor 
(spinor) fields. In such a manifestly space-time covariant presentation, one needs to choose a 
"time" direction along which to 'evolve' data specified on a 'equal time surface'. These data 
identify the configuration space variables and their velocities. The Hamiltonian formulations 
is then obtained from this Lagrangian formulation by passing through a Legendre transform. 
This identification of a time direction and a spatial slice on which data are to be specified is 
referred to as a "3 + 1 decomposition" . 

2.1 The 3 + 1 decomposition 

Let us assume that our would be space-time manifold is such as to admit a smooth function 
T : Ai — > R such that the T = constant level sets, generate a foliation. Different possible T- 
functions will generate different foliations. For this to be possible, we must have Ai ~ Ex E 3 . 

Now choose a vector field t^d^ which is transversal to the foliation i.e. every integral curve 
of the vector field intersects each of the leaves, transversally. Furthermore, locally in the 
parameter of the curve, the leaves are intersected once and only once. Normalize the vector 
field so that t^d^T = 1. This ensures that values of the T— functions can be taken as a 
"time" parameter which we denote as t. 

Fix a leaf S to and introduce coordinates, x a , a = 1,2,3 on it. Carry these along the integral 
curves of the vector fields, to the other leaves. This sets up a local coordinate system on 
Ai such that the normalized parametrization provides the coordinate t while the integral 
curves themselves are labelled by the {x a }. Note that there is no metric so Ai is not yet a 
space-time. We have only set up a coordinate system. 

Choose tensors g ab , N a , N on each of the leaves in a smooth manner ad define a space-time 
metric via the line element: 

ds 2 := - N 2 dt 2 + g ab (dx a + N a dt) (dx b + N b dt) . (2.1) 

Choosing g ab to be positive definite and N ^ ensures that the space-time metric g^ u is 
invertible. Its inverse is given by, 

g u = - N' 2 , g tb = N b N~ 2 , g ab = g ab - N~ 2 N a N b , g ac g cb = 5% . (2.2) 

We now have a space-time. The space-time metric is defined in terms of 10 independent 
functions and so there is no loss of generality. It is a convenient parametrization for reasons 
given below, but alternative parametrization are possible. 

It follows that, (i) The induced metric on the leaves is the Riemannian metric g ab . 

(ii) n M := <9 M T is normal to the leaves, since for any tangent vector X^d^, to S t , X^n^ = 
X^d^T = 0. Thanks to the normalization of t^d^, we have n M = (1, 0, 0, 0). 

(iii) n M := g^ v n v n^n^ = g tt = —N 2 < and therefore the normal is time-like and hence 
the leaves are space-like. The iVW is a unit time-like vector. 
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(iv) The original transversal vector field can be decomposed as t M = an^+N^ where N^n^ = 
and hence iV M is tangential to the leaves and N° = 0. This decomposition refers to iV 2 as the 
lapse function and iV M as the shift vector. Next, t^n^ = 1 =>■ a = —N 2 . The integral curve 
equation, d t x» = -N 2 n^ + implies for fi = a, N a = N 2 n a = N 2 g at = N 2 (N~ 2 N a ) = N a . 
This identifies the N a with the shift vector (which is spatial). 

The particular parametrization of the space-time metric can be said to be adapted to the 
pre-selected coordinate system. Since the coordinate system is defined without any reference 
to any metric, we can similarly parametrize other tensor fields, notably the co-tetrad, e 1 . 

Note: We chose an arbitrary foliation (through an arbitrary choice of a "Time" function 
and then a transversal vector field to enable us to choose coordinate on the manifold. The 
foliation provides us with a normal and the transversal vector field can be parametrized 
in terms of this normal, a lapse function and a shift vector. Varying the lapse and shift 
varies the transversal vector field relative to the foliation. If we also change the foliation, 
then the normal changes and so must the shift vector. The changes induced by lapse and 
shift correspond to making a space-time diffeomorphism and every infinitesimal space-time 
diffeomorphism can be generated by infinitesimal changes in the lapse and shift. 



2.2 Digression on tetrad formulation 

General relativity is formulated as theory consisting of tensorial fields on a manifold and a 
second rank, symmetric, non-degenerate (invertible) tensor field, g^ encoding gravitational 
phenomena. To do differential calculus on general tensor fields one also needs to define 
a covariant derivative, V M which involves the introduction of an affine connection, T x I1U , 
which is usually taken to symmetric and metric compatible i.e. V\g^ u = 0. The non-abelian 
gauge theories already introduce quantities which are not just tensors with respect to general 
coordinate transformations but also transform under the action of "an internal" group, eg 
a Higgs field <3> a , a YM potential A a , its corresponding field strength, F" etc. The index a 
indicates a response to the (adjoint) action of a group such as SU (N). Developing calculus for 
such quantities, also needs a gauge covariant derivative and a corresponding gauge connection 
eg A«. 

Consider now a quantity, e 1 (x) where \i responds to a general coordinate transformation (e 1 
transform as a covariant rank 1 tensor) and the index / responds to the local action of the 
pseudo-orthogonal group, SO(l, 3) under the defining representation. This quantity can also 
be thought of as a 4 x 4 matrix and we will take it to be an invertible matrix. This is referred 
to as a co-tetrad while its inverse quantity, is referred to as a tetrad: e^ej = 5j, e^e^ = 8 V . 
It is possible to formulate the theory of gravity in terms of a (co-)tetrad as follows. 

(1) Let e 50(1,3) i.e. A^A^ XL = r] IJ holds where 7] IJ = diag(-l, 1, 1, 1). Then the 
co-tetrad transforms as: 

r)r' u 

(e%(x'(x)) := ^{x)—ei{x) . (2.3) 

The A— transformation is referred to as a Local Lorentz Transformation (LLT) while x — > 
x'(x) is the General Coordinate Transformation (GCT). Clearly to have the derivatives of 
the co-tetrad to transform covariantly under both sets of transformations, we need two 
connections: an affine connection (not necessarily symmetric in the lower indices) and a 
Spin connection, IJ . The spin connection is anti-symmetric in the 1.1 indices and thus 
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transforms as the adjoint representation of the pseudo-orthogonal group. The derivative 
covariant with respect to the LLT is denoted by D^, that with respect to GCT is denoted 
by the usual V M while the one with respect to both will be denoted by V^. The Lorentz 
indices will be raised/lowered using the Lorentz metric i] IJ , rju, where rf K r\Kj = o~j- 



Armed with the tetrad, the spin connection and the Lorentz metric, define the following 
quantities: 



Torsion : 


T 7 (e,w) := rfe / + w / J Ae J 
T 1 av = d^ei + Wa 1 jei - (// v) 


Curvature : 


R ij {uj) := du IJ + U 1 K Au KJ 

R I J ^ = d^ v IJ + w M 1 K u u K 3 - ((j, <->■ u) 


Bianchi Identity : 
Cyclic Identity : 


{DR) 1 j = := dR 1 j + co 1 K AR K j + ooj K AR 1 K 
(DTY = R'jAe- 1 := dT 1 + oo 1 j A T J 


Metric : 


9nv{e,rj) := e^e^rju 


Christoffel Connection : 




Affine connection : 
where, 


T A ^(e, u) := | ^ | (g(e)) + \g Xa (T afMU - T, ua - T UfJi0l ) 
T aflu (e,u) := ei Q T I ^(e,w) 


These imply 

Matric compatibility : 
tetrad compatibility : 


V A ^ = 

V»4 = := V^el + u/ jei-T\ u e[ 


Compatibility =>- 


R a x^(T) = efe XJ R IJ ^(u) 



No assumption about the torsion tensor is made. 



(2) It is possible to invert the torsion equation to 'solve for' the spin connection in terms 
of the tetrad, its derivatives and the torsion tensor. All one needs to do is manipulate the 
combination Tx^u + T^v\ — T u \^ and use the invertibility of the (co-)tetrad. The result is: 

co/ J := tu/ J (e) + K/ J (e,T) (2.4) 
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w„ IJ ■= \ [e uI {0,ei - 8 v el) - e» J (d»e[ - ft,ej) - e^e XJ (ft,ef - 9 A ef ) e^] (2.5) 
K/ J := - l -e vI e XJ (T u ^ + T^ u -T^ x ) (2.6) 

The K is called the con-torsion tensor and u> is the torsion-free spin connection which is 
explicitly determined by the tetrad. The Affine connection equation is the corresponding 
inversion of the metric compatibility condition (covariant constancy of the metric) to express 
the general affine connection in terms of the torsion free Christoffel connection plus the 
torsion combinations. 

Notice that a priori, we have two connections: the affine and the spin. Both define corre- 
sponding and independent torsions (T 1 ^ and the antisymmetric part of r A fiu ). The intro- 
duction of the metric as the 'square' of the co-tetrad and the two compatibility conditions 
together identify the two torsions. 

The last equation demonstrates that we can use the tetrad and the co-tetrad to convert 
the Lorentz indices and the general tensor indices into each other with the compatibility 
conditions ensuring the two distinct curvatures also going into each other. 

Although we have referred to only 4 dimensions and Lorentz signature metric, the definitions 
generalise to any dimensions and any signature. 

(3) Four dimensions have additional features available. One can define internal dual (Lorentz 
dual) for anti-symmetric rank-2 Lorentz tensors apart from the usual Hodge dual (space-time 
dual) for 2-forms. 

Let, S IJKL and £^ va p denote the Levi-Civita symbols; These are completely antisymmetric 
in their indices and we choose the conventions: £ 0123 = 1 = £ txyz . The indices on these are 
raised and lowered by the Lorentz and the space-time metric respectively. Using these we 
define: 

X IJ := \e IJ KL X KL (Internal Dual) (2.7) 
:= \S^X^ (Hodge Dual) (2.8) 

(4) From the tetrad and the spin connection, the following Local Lorentz invariant four 
forms can be constructed whose integrals are candidate terms for an action. 

1. Hilbert- Palatini: C HP (e,oj) := \£ijklR IJ \u) A e K A e L . 

The variational equations following from this are equivalent to the Einstein equations. 
The spin-connection equation implies that the torsion vanishes and the tetrad equation 
implies the vanishing of the Ricci tensor. The Hilbert-Palatini action is thus classically 
equivalent to the Einstein-Hilbert action of the metric formulation. 

This term taken as an action with the tetrad and spin connection treated as indepen- 
dent variables is sometimes referred to as the tetrad formulation of gravity. 

2. Cosmological Constant: C\(e) := ^Sijkl^ 1 A e J A e K A e L 

This is the usual cosmological constant term, proportional to the volume form. 

3. Euler Invariant: C e (oj) := \E IJKL R 1J A R KL . 
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This 4-form is a topological term i.e. its variation under arbitrary infinitesimal changes 
in the spin connection, is an exact form and therefore the variation of its integral 
receives contributions only from the boundary values. Furthermore, explicitly, 



C E (u) = -d^-E 1 j MN co MN A (du J ! + V 



(2.9) 



4. Pontryagin Invariant: £p(u) := R IJ AR IJ . 

This 4-form is also a topological term. Furthermore, explicitly, 

C P (u) = — d |u/ j A (^du J i + k A oj k j^j | (2.10) 

The terms enclosed within the braces is the Chern-Simmons 3-form. 

5. Nieh- Yan Invariant: £ NY (e, u) := T 1 AT/ - R IJ A e 1 A e J . 

This 4-form is also a topological term which depends on both the tetrad and the spin 
connection. It vanishes if torsion is zero (for zero torsion, the second term vanishes by 
the cyclic identity.) Explicitly, 

C NY = d{e I AT 1 }. (2.11) 

Note that we have 5 different, Lorentz covariant 2-forms: T 1 , S /J := e 1 A e J , T, IJ , i? /J , i? /J . 
From these, we can form the six Lorentz invariants: T 2 ,S 2 (= = — S 2 ),S A S,i? 2 (= 
—R 2 ),R A .R, i? A E, R A E. If we are to get the Einstein equation (with a cosmological 
constant), then the T 2 and i? A E must be combined into the Nieh- Yan combination, 

(5) We will note a parametrization of the tetrad, adapted to the 3+1 decomposition, which 
leads to the corresponding metric decomposition. This can be derived from the identifica- 
tions: 

: -N 2 + g ab N a N b 

9 ab N b (2.12) 

9ab 



4 e 'aVu 



It follows, 



Co — tetrad : 

e\ := Nn T + N a Vi 
ei := Vl (free) 
9ab ■= V*V b J ruj 


Introduces nj 

n I n J r}u := — 1 , ^V^rju = 
is invertible; 


Tetrad : 

e\ := —N~ l m 

el := N^mN' + V? 


Defines V? 
n'V? := 

VfVj = 51 , VfV* = 5j + nm J 
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In this parametrization, the 16 variables in the tetrad have been traded with V^(12), N, N a (4) 
and ra J (4) variables with 4 conditions: n 2 = — l,n -V a — 0. The conditions can be viewed 
as 4 conditions on n 1 given freely chosen V£ or one condition on n 1 and 3 conditions on V/ 
given freely chosen spatial vector n l . 

Notice that we have the normalized normal, Nn^ defined by the foliation. From this we can 
construct an internal vector n 1 := ej(iVn M ). In the parametrization, we have also introduced 
an internal normalized time-like vector n 1 , determined by the freely chosen Vj . These two 
are related by the parametrization of the tetrad as, n 1 = —n 1 . 

We can view n M defined by the foliation and nj defined by a choice of Vj as two time- 
like normalized vectors in the T*(Ai). These are not identical in general and in particular 
nj is not normal to the foliation. Demanding it to be so, puts a restriction on the Vj: 
ni oc n M =>■ Hi = 0, n n° = — 1 and n ■ V a = =^ V® = 0. This implies that V/ are confined 
to T*(E). This choice is the so-called time gauge. 

Finally, we reiterate that using the tetrad and co-tetrad we can freely convert the Lorentz 
and the general coordinate indices into each other. The normalized normal (Nn^) can be 
used to define a projector, := 5% + N 2 n^n v which projects space-time tensors onto spatial 
tensors. This would lead to 3+1 decomposition (or parametrization) of all other tensorial 
quantities. 
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Chapter 3 

Symmetry Reduction 



There are different uses of the term 'symmetry reduction'. Heuristically, if S is a state space 
of a system, on which is specified an action of a group, G, which preserves the defining 
specification of the system (so that G is its symmetry group), then the space S gets "decom- 
posed" into orbits of G. The space of orbits, S/G, is 'smaller' than S and could constitute 
a simplification. S/G, is thought of as a symmetry reduction of S by G. Alternatively, one 
could restrict to the subset of the so called invariant states which may be thought of as a 
collection of trivial orbits. In our context, we will be using the term in the latter sense. The 
system could be classical or quantum mechanical. 

For example, if S is the quantum mechanical state space of a particle with a rotationally 
invariant Hamiltonian, then the subspace of the invariant states would be all the states with 
zero angular momentum. If it is the phase space of a particle with a rotationally invariant 
dynamics, then the only invariant 'state' is the origin of the configuration space with zero 
momentum. If however, S denotes the space of field configurations on a manifold, then the 
subset of invariant configurations is non-trivial. If the quantum mechanical state space of a 
system consists of distributions on a space of 'test functions', then invariant states could be 
defined as those distributions whose support consists of invariant test functions. 

If one obtains a reduction by restricting to invariant states (and invariant observables) of 
a quantum system, one has followed the first quantize, then reduce route and the reduced 
system can be thought of as a symmetric sector. This is not always possible, since one does 
not have adequate explicit control over the quantum system. Alternatively, one can consider 
invariant subspace of a classical phase space and construct a corresponding quantum theory. 
This is the first reduce, then quantize route. In general, the relation between these two 
approaches is unclear. 

While the former approach is more desirable, in practice, it is the latter approach which is 
followed commonly. We will also follow this approach. However, we will follow the methods 
- basic variables, construction of quantum Hilbert space etc - used in the full theory. The 
viability of these simplified models are then thought to constitute a test of the methods and 
premises of the full theory. The reduction of the classical theory is carried out by requiring 
certain symmetries to be exactly realized. 
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3.1 Symmetry Reduced Models 



We are already familiar with use of symmetries to simplify a problem. For example, as- 
suming spherical symmetry we choose coordinates and metric components to simplify the 
Einstein equation and obtain the Schwarzschild solution or using homogeneity and isotropy 
one obtains the FRW solutions. Thus symmetry groups (isometries) allow us to classify 
suitable ansatz for the basic variables of the theory. Note however that we are not interested 
in solving Einstein equations, but rather in obtaining a classical action with fewer degrees 
of freedom and constructing a corresponding quantum theory. In the context of spherical 
symmetry for example, this corresponds to restricting to only spherically symmetric form 
of 3-metrics: ds 2 = A 2 (t,r)dr 2 + R 2 (t,r)(d9 2 + sin 2 9d(p 2 ) and reducing the Einstein-Hilbert 
action to get an action in terms of the two field degrees of freedoms - A(r),R(r). Such 
reductions of degrees of freedom is termed mini-superspace model if the degrees of freedom 
is finite and a midi-superspace model, if the degrees of freedom is still infinite i.e. a lower di- 
mensional field theory. The former occur in homogeneous cosmologies while examples of the 
latter include spherical symmetry, certain inhomogeneous cosmological models such as the 
Gowdy models, Einstein-Rosen waves etc. Needless to say that the midi-superspace models 
are still very complicated. We will concentrate on the mini-superspace models and specifi- 
cally on (spatially) homogeneous cosmologies. We begin by defining spatially homogeneous 
space-times which are not necessarily solutions of Einstein equation. 



3.1.1 Spatially homogeneous models 

A four dimensional space-time is said to be spatially homogeneous if (a) it can be foliated 
by a 1-parameter family of space-like hypersurfaces, E t and (b) possessing a (Lie) group of 
isometries such that for each t and any two points p, q G E t there exist an isometry of the 
space-time metric which maps p to q. The isometry group G is then said to act transitively 
on each of the E t . If the group element connecting p, q is unique, the group action is said to 
be simply transitive (otherwise multiply transitive). Spatially homogeneous space-times are 
further divided into two types. 

A spatially homogeneous space-time is said to be of a Bianchi type if the group of isometries 
contains a subgroup (possibly itself), G*, which acts simply transitively on E 4 otherwise it 
is said to be of the Kantowski-Sachs type (interior of Schwarzschild solution). It turns 
out that except for the special case of E ~ S 2 x R and G = 50(3) x R, in all other cases 
one has a Bianchi type space-time. 

Transitive action implies that there must be at least three independent Killing vectors at 
each point of E t since E t is three dimensional. But there could be additional Killing vectors 
which vanish at a point. These Killing vectors generate the isotropy (or stability) subgroup, 
H of G. Since H will induce a transformation on the tangent spaces to the spatial slices, 
it must be a subgroup of SO (3) and thus dimension of G can be at most 6 and at least 3 
since the dimension of G* is always 3. All 3 dimensional Lie groups have been classified by 
Bianchi into 9 types. The classification goes along the following lines [3]. 

A Lie algebra (or connected component of a Lie group) is characterised by structure constants 
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C 1 ]K with respect to a basis Xi, satisfying the antisymmetry and Jacobi identity namely, 
[Xj,X K ] = C T j K Xi ; C I JK = C I KJ ; ^ C N IL C L JK = , I, J, K — 1,2, 3. 

(JXFC) 

Using the availability of the Levi-Civita symbols, ^/jx, £ /JA ', £123 = 1 = £ 123 , we can write 
the structure constants as, 

C'jk = £jklC li , C IJ := M IJ + S IJK A K (3.1) 

Thus, the 9 structure constants are traded for 6 M IJ (symmetric in IJ) and the 3 Ax- This 
has used only antisymmetry. The Jacobi identity implies, M IJ Aj = 0. 

Noting that the structure constants are subject to linear transformations induced by linear 
transformations, Xj — > S/Xj, on the basis of the Lie algebra, the symmetric M IJ can 
be diagonalized by orthogonal transformations and the non-zero eigenvalues can be further 
scaled to ±1: M IJ = n 1 ^ 13 . The condition M IJ Aj = implies that either Aj = (Class A) 
or Aj 7^ (class B) in which case M IJ has a zero eigenvalue and we may take the non-zero 
eigenvector Aj to be along the "1st" axis, i.e. Aj = aSi^ and n l = 0. This leads to, 

[Xj, X K \ = u'SukXj + XjA K - X K Aj . 

In the class A, there are precisely 6 possibilities organized by the rank of the matrix - 0, 1, 
2, 3 and signature for ranks 2, 3 viz (++, H — ) and (+ + +, + H — ). The eigenvalues of M IJ 
can be taken to be n 1 = ±1, 0. 

In the class B, the rank of M IJ cannot be 3 and the possibilities are restricted to the ranks 0, 
1, 2 and signatures (++, H — ) for rank 2. If the rank of M is 0, all three eigenvalues are zero 
and scaling X\, we can arrange a = 1. For rank 1, taking n% to be the non-zero eigenvalue, 
scaling Xi,X^ ensures a = 1. For rank 2 however, {712 = ±1,713 = ±1)> 110 scaling can 
preserve n 2 ,n 3 and set a = 1 (of course a = 1 is possible). 

Here is a table of the classification of Riemannian, homogeneous 3-geometries[3]: 



Type 


a 


Til 


n 2 


ns 


Remarks 


Class A 












I 














Euclidean space 












(Leads to the Kasner space-time) 


II 





1 










VIIo 





1 


1 







VIo 





1 


-1 







IX 





1 


1 


1 


S d is a special case (with isotropy) 












(Central to BKL Scenario) 


VIII 





1 


1 


-1 




Class B 












V 


1 











H 3 a special case (with isotropy) 


IV 


1 








1 




VII a 


a 





1 


1 




III 


1 





1 


-1 


sub-case of type VI a 


VI a 


a 





1 


-1 





25 



Of interests to us are the so called class A models which are characterised by the structure 
constants satisfying C 1 u = 2Aj = 0. 

When H = 5*0(3), one has homogeneity and isotropy i.e. FRW space-times. We know that 
these come in three varieties depending on the constant spatial curvature. The spatially flat 
case is of type Bianchi I while positively curved case is of type Bianchi IX. (The negatively 
curved case is in class B, type V). 

The metrics of the general Bianchi type space-times have at the most 6 degrees of freedom 
thus constituting mini-superspaces. The spatial metrics can be put in the form: ds 2 = 
gij(t)e{e J jdx l dx J ' , where e\dx % are the so called Maurer-Cartan forms on the group manifold 
G*, satisfying de 1 = — \Cj K e J A e K . When one further restricts to diagonal gu one gets the 
so-called diagonal Bianchi models. 

Remark: One should notice that restricting to a subclass of metrics amounts to introducing 
background structures from the perspective of the full theory. In the present case, these 
structures are the symmetry group and the coordinates adapted to the group action (which 
allowed the metric to be put in the specific form). This is unavoidable and constitutes a spec- 
ification of the reduced model. From the perspective of a reduced model, these structures are 
non-dynamical, analogous to the manifold structure for the full theory and therefore do not 
automatically violate background independence. Instead, the background independence now 
means that quantization procedure should not depend the metric gu which is a dynamical 
variable. 

Our basic variables however are not the 3-metric and the extrinsic curvatures. They are the 
SU(2) connection and the densitized triad. In the metric variables, the natural notion of 
symmetry is isometry while in the connection formulation it is the group of automorphisms 
of the SU(2) bundle. Thus, the cosmological models will now be understood to be charac- 
terised by groups of automorphisms of the SU(2) bundle which acts on the base manifold S 
transitively. The task is to characterise the connection and triad variables which are invari- 
ant under the group action (just as isometries mean invariant metrics). This requires more 
mathematical machinery and we will only state the conclusions^]. 

For the Bianchi models, the invariant connections and densitized triad are of the form: 

4(t,x) := SK*W(*) » E?M ■= y/m(x)pi(t)Xf(x). (3.2) 

In the above equation, a refers to spatial coordinate index, % refers to the adjoint representa- 
tion of SU(2) and / refers to the adjoint index of the Lie algebra of the symmetry (sub) group 
G* (and hence takes 3 values). The u^dx a are the Maurer-Cartan 1-forms (left-invariant 1- 
forms) on E t identified with the group manifold while Xf-£^ are the corresponding invariant 
vector fields dual to the 1-forms, i.e. a/(Xj) = co^j = <v- The go(x) is the determinant 
of the invariant metric on the symmetry group and provides the necessary density weight. 
It is regarded as a fiducial quantity and will drop out later. All the coordinate dependence 
resides in these forms, vector fields and the fiducial metric while the coefficients containing 
the t dependence are the basic dynamical variables^. 

lr The remaining, Class B models are thought not have a Hamiltonian formulation and hence are not 
amenable to analysis by canonical methods [3]. 

2 A few essential details from Forgacs and Manton are summarised in the appendix. 

3 Similar decomposition is made for all quantities with spatial and the Lorentz indices. The contravariant 
spatial index is expressed using the invariant vector fields and the covariant one using the invariant 1-forms. 
In particular, K l a := K\w{ , T l a := 
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If we have isotropy in addition, then the degrees of freedom are further reduced: <&} := 
cA}, p\ := pA{ and there is only one degree of freedom left. Here the A's are a set of 
orthonormal vectors satisfying, A}Aj = 5jj, A l jAjA K €ijk = zuk- The phase space variables 
c,p are gauge invariant. 

The intermediate case of diagonal models arises from a choice := cjA),pl := p A\ (no 
sum over I). The residual (SU(2)) gauge transformations act on the A's and leaving the 
ci,p ! as the gauge invariant phase space variables thereby solving the Gauss constraint at 
the outset@. Thus there are only 3 degrees of freedom [5]. 

Having identified relevant degrees of freedom parameterising quantities invariant under sym- 
metry transformation, the next task is to obtain the symplectic structure (basic Poisson 
brackets) and simplify the expressions for the constraints. 

Symplectic form: In the full theory, this is given by (87rG7) _1 J s d 3 xA l a (t, x)E®(t, x). Direct 
substitution gives, 

— / d i xAl(t,x)E^t,x) = —¥ lP l If d s x^\ , =► {*j,j£} = ^5}#. (3.3) 
«7 is «7 Us J Vq 

The quantity in the braces is the fiducial volume, Vb, of E t . For spatially flat, isotropic case, 
the slice is non-compact and the fiducial volume is infinite. This problem is addressed by 
restricting to a finite cell whose fiducial volume is finite. One has to ensure that the final 
results do not depend on the fiducial celJE The dependence on the fiducial volume is gotten 

]/3 2/3 

rid off by redefining the basic variables as $ — > §V , p — >■ pV . If we have isotropy, 
the symplectic form would become -^Vodp which leads to (after rescaling) to the Poisson 
bracket, {c, p} = With this rescaling understood, we will now effectively put Vq = 1. 

Curvature: The curvature corresponding to the invariant connection above, is obtained as: 

F l := dA l + -e l lk A 3 A A k := -F) k uj j A u K (3.4) 
2 2 

.-. F) K = -QiC 1 JK + e* jk &j& K (3.5) 



Gauss Constraint: The full theory expression is: 
G(A) , 



A 2 



Pi 



-V C J U 



Vn 



Gi = (K7) _1 {-PIC J U + ej&jpl} ■ 



(3.6) 
(3.7) 



Notice that for the class A models, the first term is zero and for the diagonal models the 
second term vanishes as well (since e is antisymmetric in j, k while the A factors are symmetric 
in j, k). There are no continuous gauge invariances left. Note that the first term in eqn (13. 7p . 

4 There is still a discrete invariance remaining and involves changing the sign of two of the triad and 
connection components. 

5 This is discussed in more details in section 14.2.21 
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is a surface term which could vanish if E has no boundaries. But this would not be true for 
say spatially flat models which will have the £ as a cell on which the invariant vector fields 
need not vanish. The integrand however is proportional to C J 13 and these vanish for the 
class A models. 

Diffeo Constraint: 

C diS (N) := — / N a (x)E^x)F^ ab (x) - [ N a (x)Al l (x)G t (x) 

N a {t,x) := N'tyXftx) (3.8) 

Afl 

:. N J Cj = — [{C K j K ¥ I + C K IJ ^ K )p(\ (3.9) 

This constraint again vanishes for diagonal, class A models. 

Hamiltonian Constraint: The full theory Hamiltonian constraint is given by, 

CnUN) ■= y k I N ^^\ - 2(1 + ^ K l K i\\ ( 3 - 10 ) 

To carry out the integration, we need to note the expressions: 

VFo = \.eij K t ahc ^l^ , -U = e abc e IJK X?X b jX c K . (3.11) 



This leads to, 

detg = det(Et) := ^ ahc ^\g^l 2 X a 1 X)X^p\p]pK = i^W^ (3-12) 
Now the integration can be carried out immediately to give, 

N 



#grav ^ 



P * P ' J> kF k j j - 2(1 + 7 2 )^M ] } 



(3.13) 



In the above, = 7~ 1 ($} — T}). These expressions are valid for general Bianchi models. 

At this stage, we could in principle attempt to carry out the usual Schrodinger quantization 
with being multiplicative operators and P/ being the derivative operators. Both transform 
covariantly under the action of SU(2). 

However, we can also imagine 'specializing the holonomy-flux variables' of the full theory, 
for these symmetric fields. It is natural to choose edges along the symmetry directions 
i.e. along integral curves of the Xf vector fields. It follows that due to homogeneity, the 
path ordered exponentials, holonomies, become just the ordinary exponentials, hi($) := 
h ei (&) := "Pexp{ f Q'jTiUia dx a } = exp{$}(t)rj j^^^dx 0- }. There is no sum over I in these 
expressions. These can be further expressed using the identity e %e ' n ' (T = cos(0)+ in -asm9. The 
holonomy is then given in terms of 9 ~ ^JW^\ which is gauge invariant and two angular, 
gauge variant components corresponding to the direction n ~ unit vector in the direction 
of $}. A simplification occurs if we further restrict to the diagonal models: := C/A} 

— * 

which makes the n = Aj and now the matrix elements of these holonomies can be obtained 
from the elementary functions, e^ I ~>° 1 ^ 2 . These have been termed as the point holonomies. 
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The fluxes through surfaces perpendicular to the symmetry directions, likewise simplify to 
E SjK (f) = p J Af/ 1 J s 'go(x)e a b c XfdS bc oc p 1 . Unlike the flux operators in the full theory, 
these fluxes Poisson commute among themselves. Thus, in the diagonal models, we can 
extract gauge invariant phase space coordinates, with the holonomies and fluxes having the 
usual Poisson algebra. In quantum theory, a useful triad representation can then be set-up. 

Point holonomies and commuting flux variables are new features which arise in the (diagonal) 
mini-superspace reduction. These are also responsible for the relative ease of analysis possible 
for these models. This will be discussed more below. 

What about inhomogeneous models? There have fewer efforts regarding these. Among the 
inhomogeneous models, the reduction for the Gowdy model on 3-torus can be seen in [6], 
while spherical symmetric model can be seen in [TJ. Martin's lattice model is briefly discussed 
in the appendix. 
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Chapter 4 

Singularity Resolution in Quantum 
Theory 

The most detailed analysis of the singularity resolution is available for the homogeneous and 
isotropic geometry coupled to a massless scalar and this is the case that we discuss below. 
Prior to 2005, the singularity resolution was understood as the deterministic nature of fun- 
damental equation (the Hamiltonian constraint) and in terms of an effective picture deduced 
either from the WKB approximation of by taking expectation values of the Hamiltonian. 
In this sense, resolution of singularities was seen for (i) FRW coupled to a scalar field with 
arbitrary positive semidefinite potential and (ii) diagonalised Bianchi class A (anisotropic) 
models. These resolutions were seen as an implication of the inverse triad quantum correc- 
tions which were present in the matter sector (and in the curvature for non-flat models). 
Post 2005, it was realized, at least for the FRW case, that the holonomy corrections by 
themselves could also resolve singularities. This required restriction to massless scalar and 
treating it as a clock variable, thereby paving the way for construction of physical states, 
Dirac observables and physical expectation values. Although restricted to special matter, it 
allows completion of the quantization program to the physical level and throws light on how 
a quantum singularity resolution may be viewed. For this reason, we this case is discussed 
in detail. Subsequently, Madhavan also showed another quantization for the same case, also 
completed to physical level, wherein holonomy corrections are absent and singularity reso- 
lution is achieved by inverse triad corrections only. There are also some issues which have 
been better understood in the past few years. These are briefly summarised and discussed 
in sections KT2\ and l4T2~31 . 



4.1 FRW, Classical Theory 

Classical model: Using coordinates adapted to the spatially homogeneous slicing of the 
space-time, the metric and the extrinsic curvature are given by, 

ds 2 := -dt 2 + a 2 (t) {dr 2 + r 2 dtt 2 } := - dt 2 + a 2 (t)ds 2 comoviQg . (4.1) 
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Starting from the usual Einstein-Hilbert action and scalar matter for definiteness, one can 
get to the Hamiltonian as, 
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Thus, H — -<->■ Friedmann Equation. For the spatially flat model, one has to choose a 
fiducial cell whose fiducial volume is denoted by Vq. 

In the connection formulation, instead of the metric one uses the densitized triad i.e. instead 
of the scale factor a one has p, \p\ := a 2 while the connection variable is related to the 
extrinsic curvature as: c := 7a (the spin connections is absent for the flat model). Their 
Poisson bracket is given by {c, p} = (87^7) / (3 Vo). The arbitrary fiducial volume can be 
absorbed away by defining c := Vq 1//3 c, p := V 2 ^ 3 p. Here, 7 is the Barbero-Immirzi parameter 
which is dimensionless and is determined from the Black hole entropy computations to be 
approximately 0.23 [8]. From now on we put 8ttG := n. The classical Hamiltonian is then 
given by, 
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For future comparison, we now take the potential for the scalar field, V(4>) to be zero as well. 

One can obtain the Hamilton's equations of motion and solve them easily. On the constrained 
surface (H = 0), eliminating c in favour of p and p^, one has, 
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(4.6) 
(4.7) 



Since is a monotonic function of the synchronous time t, it can be taken as a new "time" 
variable. The solution is determined by p(4>) which is (i) independent of the constant p^ and 
(ii) passes through p = as — > ±00 (expanding/contracting solutions). It is immediate 
that, along these curves, p((p), the energy density (p~ 6 p^/2) and the extrinsic curvature 
diverge as p — )■ 0. Furthermore, the divergence of the density implies that d>(t) is incomplete 
i.e. t ranges over a semi-infinite interval as <p ranges over the full real linqj. 



1 For the FRW metric, integral curves of dt are time-like geodesies and hence incompleteness with respect 
to t is synonymous with geodesic incompleteness. 
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Thus a singularity is signalled by every solution p(4>) passing through p = in finite syn- 
chronous time (or equivalently by the density diverging somewhere along any solution). A 
natural way to ensure that all solutions are non-singular is to ensure that either of the two 
terms in the Hamiltonian constraint is bounded. Question is: If and how does a quantum 
theory replace the Big Bang singularity by something non-singular? 

There are at least two ways to explore this question. One can imagine computing corrections 
to the Hamiltonian constraint such that individual terms in the effective constraint are 
bounded. This approach presupposes the classical framework and thus will have a domain 
of validity of these corrections. Alternatively and more satisfactorily, one should be able to 
define suitable Dirac observables whose expectation values will generate the analogue of p(<p) 
curves along which physical quantities such as energy density, remain bounded. Both are 
discussed below. 



4.2 FRW, Quantum Theory 

Schrodinger Quantization: In the standard Schrodinger quantization, one can introduce 
wave functions of p, <f> and quantize the Hamiltonian operator by c — > itiK'j/Zdp , p^ — > —ihd^, 
in equation (I4.5p . With a choice of operator ordering, H^/(p, 0) = leads to the Wheeler-De 
Witt partial differential equation which has singular coefficients. We will return to this later. 

Loop Quantization: The background independent quantization of Loop Quantum Gravity 
however suggest a different quantization of the isotropic model. One should look for a 
Hilbert space on which only exponentials of c (holonomies of the connection) are well defined 
operators and not c. Such a Hilbert space is obtained as the representation space of the 
C* algebra of holonomies. In the present context this algebra is the algebra of almost 
periodic functions of c, finite linear combinations of functions of the form e lXjC , Xj G K. 
Inner product (analogue of the Ashtekar-Lewandowski measure) on the space of the almost 
periodic functions is given by: 

r-T 



:= lim -L f dc¥*( C )*(c) . (U 

1 — >oo ZJ. J —T 



The single exponentials form an orthonormal set. Let us denote it as, (c\fj) :=exp{|/ic}, \i G 
M. The holonomy-flux representation can now be made explicit as: 

o 

h v \fi) := e^ c |/x) = + (4.9) 

Notice that that the triad operator has every real number as a proper eigenvalue (i.e. has 
a corresponding normalizable eigenvector, the spectrum is discrete). This implies that the 
holonomy operator, is not weakly continuous in the label v i.e. arbitrary matrix elements of 
h v are not continuous functions of v. Therefore one cannot define a c operator. Note that 
the volume operator^ is given by V := |p| 3 / 2 . 

Inverse Triad Operator: The fact that spectrum of the triad operator is discrete, has a 
major implication: inverses of positive powers of triad operators do not exista These have 

2 The domain of the inverse power operator(s) will have to exclude the subspace corresponding to the zero 
eigenvalue of the triad operator. But this makes the domain non-dense and its adjoint cannot be defined [9]. 
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to be defined by using alternative classical expressions and promoting them to quantum 
operators. This can be done with at least one parameter worth of freedom, eg. 
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(4.10) 



Only positive powers of \p\ appear now. However, this still cannot be used for quantization 



since there is no c operator. One must use holonomies: hj(c) 



are 



anti-hermitian generators of SU(2) in the j th representation satisfying Tr, (rjTj) = — + 
1) (2j + l)Sij, A 1 is a unit vector specifying a direction in the Lie algebra of SU(2) and /iq 

1/3 

is the coordinate length of the edge used in defining the holonomy. It is a fraction ofV . 
Using the holonomies, 
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which can be promoted to an operator. Two parameters, /io G i and j G N/2, have crept in 
and we have a three parameter family of inverse triad operators. The definitions are: 
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All these operators obviously commute with p and their eigenvalues are bounded above. This 
implies that the matter densities (and also intrinsic curvatures for more general homogeneous 
models), remain bounded over the classically singular region. Most of the phenomenological 
novelties are consequences of this particular feature predominantly anchored in the matter 
sector. In the effective Hamiltonian computations, this modification will imply the second 
term in the Hamiltonian constraint ( 14. 5 p is rendered bounded implying singularity avoidance. 

We have also introduced two scales: p := |/xo^p an d 2jp := |/i (2j)£p. The regime 
\p\ <C po is termed the deep quantum regime, p ^> 2jp is termed the classical regime and 
Po |p| ^ 2jp is termed the semi-classical regime. The modifications due to the inverse 
triad defined above are strong in the semi-classical and the deep quantum regimes. For 
j = 1/2 the semi-classical regime is absent. Note that such scales are not available for the 
Schrodinger quantization. 

The Gravitational Constraint: Since c operator does not exist, the gravitational Hamil- 
tonian (the first bracket in eq. (l4.5p ). has to be expressed in an equivalent form using 
holonomies. For this, let us go back to the full theory Hamiltonian: 
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Now use the two identities: 



s ijk E ai E b J _ ^ 4 sgn p , . / (/lo) f / (mo) \-i 



vldetBj V «7/ioK) 1/3 



r « ,VU (4.15) 



F* = -2 lim Trl -B* ^r 1 ^ (4.16) 



A n ->0 



^(Mo) ;= ^o)^o)^(Mo)j 1 (4.17) 

In the above, the fiducial cell is thought to have been sub-divided into smaller cells of side 
I := Vq 1 ^ 3 . The area of a plaquette is Au = ^Vo- The plaquette is to be shrunk such that 
its area goes to zero. The superscript (/io) on the holonomies is to remind that the length 
of the edge is /i - The 1- forms co l a are the fiducial 1-forms whose square gives the fiducial 
metric and the e abc is the (fiducial) metric dependent Levi-Civita density. 

In quantum geometry however there is a gap in the spectrum of area operator and thus it is 
not appropriate to take the area to zero, but at the most to the smallest possible eigenvalue. 
Independently, if we force the limit, it will imply /io — > which in turn amounts to defining 
c operator which does not exist on the Hilbert space. 

Substituting these in the (14.141) and carrying out the integration over the cell leads to (sup- 
pressing the /io superscript on the holonomies), 

tfgrav = ~ ^ 3 E^ Tr (MA r VWW) (4-18) 

In the above, we have used j = 1/2 representation for the holonomies and V denotes the 
volume function. In the limit /x — > one gets back the classical expression Q 

If we promote this expression to a quantum operator (modulo ordering ambiguities) on the 
LQC Hilbert space constructed above, then we cannot take the limit /i — > because it would 
imply that c exist which we have shown to be impossible. Thus, at the quantum level we 
should not take the limit /i — > 0. The best we can do is to take reduce /io such that the area 
reaches its smallest possible (and non-zero due to the gap) eigenvalue A := (2\/37r7)£p. But 
which area do we consider, the fiducial or the physical? These are related by a factor of \p\. 
It seems appropriate to choose the physical area, which implies that we must take /io to be a 
function Jx of p given by, p(p) := ^ A/\p\. Note that this is one prescription to interpret the 
limitation on shrinking of the plaquette. There are others which will be mentioned later. In 
the following we will continue to use the /i notation and replace it by p,(p) when needed. 

While promoting this expression to operators, there is a choice of factor ordering involved 
and many are possible. We will present two choices of ordering: the non-symmetric one 
which keeps the holonomies on the left as used in the existing choice for the full theory, and 
the particular symmetric one used in [10J. 

= -^-W (sin^cos^f - costfv sin^) (4.19) 



3 The expression for the Hamiltonian constraint follows exactly from the full theory procedure. Starting 
from the equation (|4.14|) . the integral will be replaced by a sum over smaller cells of a triangulation. The 
size parameter of the cells will drop out thanks to density weight 1 of the Hamiltonian. Due to homogeneity, 
contribution from each cell will be the same and hence the total sum will be number of cells of the triangula- 
tion times the contribution of one cell. There are exactly /Xq 3 cubical (say) cells with side of length /ioT^ 3 
and this produces the factor of /i^~ 3 in equation (|4.18j) . The Vb of course disappears as in the full theory 
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/Vsym 24«(sgll(p)) . / . (iQC y /i C HQC y . /X C\ . 

H lrZ = ^3^2 sm/xoc (sin— Vcos— - cos— Vsm— j sm/x c (4.20) 

At the quantum level, /xo cannot be taken to zero since c operator does not exist. The action 
of the Hamiltonian operators on |/x) is obtained as, 

tf™r yr >> = ^^(W-^o)(l/^ + 4/io)-2|/i) + |/x-4/io» (4.21) 
3 

#gravl^> = ^3g| [1^+3^0 ~ ^Wo I I A* + 4 /i > + I V» ~ ^-3 M I I A* — 4 ^0> 

- {|U - Wl + IV W - ^-3 M |} (4-22) 
where V M := (|7^p|/i|) 3//2 denotes the eigenvalue of V. 

We also have the Hilbert space for the matter degrees which for us is a single scalar, and 
the full kinematical Hilbert space is the tensor product of the L 2 (MBohr, c^Bohr) ® ^matter- 

Wheeler-DeWitt Difference Equation: Let |\&) := ^(a*, 0)1/-*), where the sum is 
over a countable subset of R, the coefficients \I/ (/x, 0) are valued in the matter Hilbert space 
and the argument is a reminder of that. The Hamiltonian constraint is imposed on these 
|^} which leads to the Wheeler-DeWitt equation for the coefficients. Thanks to the presence 
of the trigonometric operators, this equation is a difference equation. In the Schrodinger 
quantization, this would be a differential equation. 

For the non-symmetric operator we get, 

A(fi + 4/i )^(/i + 4/i , 0) - 2A(/x)*(/x, 0) + A(n - 4/i )^(/i - 4/i , 0) 

= - yA*o7 3 4^ma«er-(A*)*(M,0) (4.23) 

where, A(fj) := V^ +W) — V^_ Mo and vanishes for /x = 0. 
For the symmetric operator one gets, 

f (fi + 4/x , 0) + fo(n)*(», <P) + - 4^o, 0) 

= - —i2 3 ol 3 e 2 P H matter (fi)^(i2,(j)) where, (4.24) 

: = 11^+3^,-1^+^1 , /_(/x) := /+(/x-4/x ), / := - f+(fi) - f-(fi) . 

Notice that /+(— 2/i ) = = /_(2/x ), but /o(a0 is never zero. The absolute values have 
entered due to the sgn(p) factor. 

The difference equations relate ^(aO's only for /i's in a "lattice" , £^ := {/i = /x+4/x n, n G Z} 
and the coefficients labelled by different lattices are completely independent. The /x G 
[0,4/i ), label different superselected sectors. 

The equations are effectively second order difference equations and the \I/(/x,0) are deter- 
mined by specifying ^ for two consecutive values of /x eg for /x — (i + 4/i iV and /x = 
/x + 4/xo(iV + 1). Since the highest (lowest) order coefficients vanishes for some /x, then the 
corresponding component \l/(/x,0) is undetermined by the equation. Potentially this could 
introduce an arbitrariness in extending the \l/ specified by data in the classical regime (eg 
/x ^> 2j) to the negative /x. Potentially, maintaining determinism of the quantum wave 
function, is one of the restrictive criteria for choosing the ordering. 



36 



For the non-symmetric case, the highest (lowest) A coefficients vanish for their argument 
equal to zero thus leaving the corresponding \1/ component undetermined. However, this 
undetermined component is decoupled from the others. Thus apart from admitting the trivial 
solution ^(jti, <fi) := $(0)5^0, V/i, all other non-trivial solutions are completely determined 
by giving two consecutive components: ^(fi, 4>),^(fi + 4yU , <f>)- 

For the symmetric case, due to these properties of the /±,o (/•*), it looks as if the difference 
equation is non- deterministic if \i = 2/i + 4/i ^, n G Z. This is because for // = — 2/i , 
^>(2fi ,(j)) is undetermined by the lower order \l/'s and this coefficient enters in the deter- 
mination of \I/(2/!o,0). However, the symmetric operator also commutes with the parity 
operator: (Tl\I/)(yU, (f>) := \&(— /i, (f>). Consequently, ^/(2fi ,(f)) is determined by \&(— 2/z , 0). 
Thus, we can restrict to /i = 2/i + 4fc/i , > where the equation is deterministic. 

In both cases then, the space of solutions of the constraint equation, is completely determined 
by giving appropriate data for large \fi\ i.e. in the classical regime. Such a deterministic 
nature of the constraint equation has been taken as a necessary condition for non-singularity 
at the quantum level @. 

Effective Hamiltonian: By introducing an interpolating, slowly varying smooth function, 
^(p(fi) := ^jipfi), and keeping only the first non- vanishing terms, one deduces the Wheeler- 
De Witt differential equation (with a modified matter Hamiltonian) from the above difference 
equation. Making a WKB approximation, one infers an effective Hamiltonian which matches 
with the classical Hamiltonian for large volume (/i 3> /^o) and small extrinsic curvature 
(derivative of the WKB phase is small). There are terms of o(h°) which contain arbitrary 
powers of the first derivative of the phase which can all be summed up. The resulting 
effective Hamiltonian now contains modifications of the classical gravitational Hamiltonian, 
apart from the modifications in the matter Hamiltonian due to the inverse powers of the 
triad. The largest possible domain of validity of effective Hamiltonian so deduced must have 

\p\ >p mm- 

An effective Hamiltonian can alternatively obtained by computing expectation values of the 
Hamiltonian operator in semi-classical states peaked in classical regimes [13] . The leading 
order effective Hamiltonian that one obtains is (spatially flat case): 



H. 



non— sym 
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ZM/;)snr(/,,,r)- ( A(p) - ^B+(p) 



matter j 



eff 16ttG* V/^o7 3 4 

B+(p) := A(p + Ap ) + A(p-A Po ) , A{p) := {\p + p \ 3/2 - \p - Po\ 3/2 ) , (4.25) 

X 1 

p := -74/* , Po ■■= «74a*o ■ 
o o 

For the symmetric operator, the effective Hamiltonian is the same as above except that 
B+(p) -> f+(p) + f-(p) and 2A(p) ^ f + (p) + f_(p). 

The second bracket in the square bracket, is the quantum geometry potential which is neg- 
ative and higher order in £p but is important in the small volume regime and plays a role in 
the genericness of bounce deduced from the effective Hamiltonian [14]. This term is absent in 
the effective Hamiltonian deduced from the symmetric constraint. The matter Hamiltonian 



4 For contrast, if one just symmetrizes the non-symmetric operator (without the sgn factor), one gets a 
difference equation which is non- deterministic. Note that this issue arises only in one superselection sector 
so may not really be an issue. However, requiring deterministic equation in all sectors could be invoked as 
a criterion to discriminate between different factor ordering. 
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will typically have the eigenvalues of powers of inverse triad operator which depend on the 
ambiguity parameters j, I. 



We already see that the quantum modifications are such that both the matter (due to inverse 
volume corrections) and the gravitational part (due to holonomy corrections) in the effective 
Hamiltonian, are rendered bounded and effective dynamics must be non-singular. 

For large values of the triad, p ^> p , B + {p) ~ Qpo^/p — o(p~ 3 / 2 ) while A(p) ~ Spo^/p — 
o(p~ 3 / 2 ). In this regime, the effective Hamiltonians deduced from both symmetric and non- 
symmetric ordering are the sam^f). The classical Hamiltonian is obtained for p — >■ 0. From 
this, one can obtain the equations of motion and by computing the left hand side of the 
Friedmann equation, infer the effective energy density. For p ^> po one obtains 0, 

3 f a 2 \ ( ggatteA 87rG^ 2 ( ggatteA ) 2 /, (A OR\ 

s^Uj :=Peff= \ ^VWJ! ' p:=a/ ( } 



The effective density is quadratic in the classical density, p c \ := -ff m atterP~ 3 ^ 2 - This modifi- 
cation is due to the quantum correction in the gravitational Hamiltonian (due to the sin 2 
feature). This is over and above the corrections hidden in the matter Hamiltonian (due to 
the "inverse volume" modifications). As noted before, we have two scales: po controlled by 
Ho in the gravitational part and 2p j in the matter part. For large j it is possible that we 
can have po <C p <C 2p j in which case the above expressions will hold with j dependent 
corrections in the matter Hamiltonian. In this semi-classical regime, the corrections from 
sin 2 term are smaller in comparison to those from inverse volume. If p ^> 2poj then the 
matter Hamiltonian is also the classical expression. For j = 1/2, there is only the p 3> po 
regime and p c \ is genuinely the classical density. 

To summarize: 

(1) The connection formulation, in the homogeneous and isotropic context, uses variables 
(c, p G R) in terms of which the classical singularity (p = 0) is in the interior of the phase 
space. By contrast, in the ADM variables (a > 0,K), in the same context, the classical 
singularity (a = 0) is on the boundary. This requires a boundary condition on the quantum 
wave functions to be specified in the deep quantum region where the classical framework is 
suspect. When the singularity is in the interior, only a continuation of the quantum wave 
function is required, given its specification in the semi-classical region. 

(2) The connection variables also strongly motivate the very different loop quantization. Its 
immediate implications are two types of corrections - the holonomy corrections and the 
inverse triad corrections. Either of these is sufficient to indicate a bounce in the effective 
Hamiltonian picture. The same use of holonomies make the Wheeler-DeWitt equation, a 
difference equation. 

(3) The analysis at the level of effective Hamiltonian already indicates (i) replacement of 
big bang by big bounce; (ii) natural prediction of an inflationary (accelerated) phase; (iii) 
singularity resolution for more general homogeneous models with curvature. 

5 The effective Hamiltonian then reduces to — - 7 _2 v /p [/i^~ 2 sin 2 (/Ltoc)] . This is also the Hamiltonian in eq. 
(|4.18[) for non-zero {1q. 

6 For p in the semi-classical regime, one should include the contribution of the quantum geometry potential 
present in the non-symmetric ordering, especially for examining the bounce possibility |12) . 
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(4) There are at least three distinct ambiguity parameters: llq related to the fiducial length 
of the loop used in writing the holonomies; j entering in the choice of SU(2) representation 
which is chosen to be 1/2 in the gravitational sector and some large value in the matter 
sector; I entering in writing the inverse powers in terms of Poisson brackets. The first one 
was thought to be determined by the area gap from the full theory. The j = 1/2 in the 
gravitational Hamiltonian seems needed to avoid high order difference equation and larger 
j values are hinted to be problematic in the study of a three dimensional model [15] . Given 
this, the choice of a high value of j in the matter Hamiltonian seems unnatural Nevertheless 
the higher values of j in the matter sector allow for a larger semi-classical regime. The I 
does not play as significant a role. 

(5) LQC being a constrained theory, it would be more appropriate if singularity resolution 
is formulated and demonstrated in terms of physical expectation values of physical (Dirac) 
operators i.e. in terms of "gauge invariant quantities". This can be done at present with 
self-adjoint constraint i.e. a symmetric ordering and for free, massless scalar matter. 

Physical quantities and Singularity Resolution: When the Hamiltonian is a constraint, 
at the classical level itself, the notion of dynamics in terms of the 'time translations' generated 
by the Hamiltonian is devoid of any physical meaning. Furthermore, at the quantum level 
when one attempts to impose the constraint as H\^f) = 0, typically one finds that there 
are no solutions in the Hilbert space on which H is defined - the solutions are generically 
distributional. One then has to consider the space of all distributional solutions, define a 
new physical inner product to turn it into a Hilbert space (the physical Hilbert space), define 
operators on the space of solutions (which must thus act invariantly) which are self-adjoint 
(physical operators) and compute expectation values, uncertainties etc of these operators to 
make physical predictions. Clearly, the space of solutions depends on the quantization of 
the constraint and there is an arbitrariness in the choice of physical inner product. This is 
usually chosen so that a complete set of Dirac observables (as deduced from the classical 
theory) are self-adjoint. This is greatly simplified if the constraint has a separable form 
with respect to some degree of freedom^]. For LQC (and also for the Wheeler-De Witt 
quantum cosmology), such a simplification is available for a free, massless scalar matter: 
-f^matter(0, P<j>) '■= Ip^,\p\~ 3 ^ 2 ■ Let us sketch the steps schematically, focusing on the spatially 
flat model for simplicity [TOf [T7]. 

1. Fundamental constraint equation: 
The classical constraint equations is: 

-^ c2 V¥\+8irGpl\p\- 3/2 = = C grav + C matter ; (4.27) 

The corresponding quantum equation for the wave function, \I/(p, </>) is: 

87rG$tf(p,0) = [Bip^C^ip^) , [B{p)\ is eigenvalue of |5p7* ; (4.28) 

Putting ^0 = —ihdff), p := ^-/i and B(p) := (^)~ 3 / 2 £>(/i), the equation can be written 
in a separated form as |§, 

:= -9(^*0,0). (4.29) 



7 For an alternative view on using large values of j , see reference [16] . 
8 A general abstract procedure using group averaging is also available. 

9 Our primary goal here is to compare the classical geometry (Wheeler-DeWitt quantization) and quantum 
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The operator for different quantizations is different. For Schrodinger quantization 
(Wheeler-De Witt), with a particular factor ordering suggested by the continuum limit 
of the difference equation, the operator 0(/x) is given by, 

e Sch (/i)*(/i,0) = -^M 3 ^0 M (V£W^)) (4.30) 

while for LQC, with symmetric ordering, it is given by, 

e LQC (/i)*(/i,0) = -[ J B(/i)]- 1 {C + (/i)vI/( /U + 4 / io,0) + C o (/i)v[/( / i,0)+ 

C + M ■= ^| |/x + 3/i | 3/2 -|/i + /io| 3/2 | , (4.31) 
C-(fi) := C + (/i-4 yUo ) , C ( M ) := -C+M-C-Qm). 

Note that in the Schrodinger quantization, the B Sch (jj) = |/i| -3 ^ 2 diverges at fx = 
while in LQC, 5lqc(a0 vanishes for all allowed choices of ambiguity parameters. In 
both cases, B(n) ~ |/i| -3 ^ 2 as |/i| — > oo. 

Inner product and General solution: 

The operator turns out to be a self-adjoint, positive definite operator on the space of 
functions \l/(/i, 0) for each fixed with an inner product scaled by B(fi). That is, for 
the Schrodinger quantization, it is an operator on L 2 (M, Bs c ^fi)d^i) while for LQC it is 
an operator on L 2 CR Bohr , B Bohr (fi)dfi Bohr ) . Because of this, the operator has a complete 
set of eigenvector^: 

0e fc (/i) = u 2 (k)eM , (e fe |e fe /) = 6(k,kf) , k, k' e R. (4.32) 

Consequently, the general solution of the fundamental constraint equation can be ex- 
pressed as 

J dk ^ + (A;)e fe ( / u)e^ ( ^ + ^_(A;)e fc (/i)e- iaj(fc) * . (4.33) 

The orthonormality relations among the ek(fj) are in the corresponding Hilbert spaces. 
Different quantizations differ in the form of the eigenfunctions, possibly the spectrum 
itself and of course u{k). In general, these solutions are not normalizable in L 2 (MBohr x 
R, -E>Bohr(/-0^/-tBohr x dfj,) , i.e. these are distributional. 

Remarks: 

(1) The operator acts in each of the superselected sector (thanks to the difference 
equation structure) and these are separable. Hence Dirac-5 appears in general when 
the label takes continuous values. 

(2) The group averaging can be seen as follows. Given any /(//, <p) one defines its group 
average, 

r°° B 2 
*fM) ■= J d\e iXC ™ f(ji,<f>) ; Oat :=— + ©, 



geometry (loop quantization). Consequently, the gravitational constraint is quantized in two different ways 
but for simplicity, the matter sector is quantized in the usual Schrodinger way. In both quantizations, 
= —ihdcf, and there is no cf> dependence in the matter Hamiltonian, so the two quantum Hamiltonian will 
have identical expressions. However, using cf> as labelling an 'emergent time' would be questionable. 

1/4 

10 For the Schrodinger quantization, the explicit eigenfunctions are: efc(/z) := iHL— e tkln \^\ an( j the eigen- 
values are: u 2 {k) := ?f{k 2 + 1/16) [17]. 
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The eigenfunctions of the constraint operator are of a product form thanks to separa- 
bility, 

—^(0) = -<tV(0) , Ge k (fi) = u 2 (k)e k (fi) , g a {<j>) = e %a \ k,a G R 



Expanding the general function /(//, 0) in the eigenbasis of the constraint operator, 
f{n,4>) := j dk J da A K(J g a (<f>)e k (fj) , 

implies, 



*f(ji,<p) = I dk j da A Ka g a (<f))e k (ii)5(a 2 - u\k)) 

dk 



oo 



2\u(k)\ 



[{A k , u{k) e^ + 4,-, W e- luW } eM] (4.34) 
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In the second equation above, we have carried out the integration over a using S(a 
oj 2 {k)) = 2 \J(k)\ (^(°" ~~ U W) + + w(k)). Clearly the group average of a general 
function reproduces precisely the general solution given in equation (I4.33p . 

Choice of Dirac observables: 

Since the classical kinematical phase space is 4 dimensional and we have a single 
first class constraint, the phase space of physical states (reduced phase space) is two 
dimensional and we need two functions to coordinatize this space. We should thus look 
for two (classical) Dirac observables: functions on the kinematical phase space whose 
Poisson bracket with the Hamiltonian constraint vanishes on the constraint surface. 
Specific values of these functions serve to label the physics states. Thus physical 
observables are values of the Dirac observables. Classically, the Dirac observables can 
be obtained as follows. 

Our constraint is: — ^ r y~ 2 c 2 p 2 ~ 0. A Dirac observable is a function f((f>,p^,c,p) 
whose Poisson bracket with the constraint vanishes on the constraint surface. We 
can describe the constraint surface by solving for c as: c(0,p^,p) := ±a/ n/Q^p^p -1 
(say). Consequently it should suffice to consider the Dirac observables to be a func- 
tion of (<fr,p$,p) (we need only two independent Dirac observables). Then the Dirac 
observables are defined by the differential equation, 

Evidently, / = p^ is a Dirac observable. For the second one, we can choose any 
function of (. A particularly convenient choice is: f(C(4>,p)) '■= #exp{=FA/2/?/3C} = 
p exp{=|=A/2fi;/3(</> — 0o)}- These Dirac observables taking a particular value, say p*, 
define curves in the (p, (j)) plane, p(4>) := p*exp{±A/2ft/3(0— 4> )} which are the classical 
solutions in (14. 7p . 

In the quantum theory, the notion of Dirac observable is that it is an operator which 
maps solutions of the constraint to (other) solutions. We already have the general 
solution in eq f 14 . 3 3 1) which is obtained via unitary evolution (in 0) from an ini- 
tial ty±(4>o, /•*)• Hence a Dirac observable is constructed by defining an operator on 
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\l/±(0o,/i) to generate a new 'initial state', and evolving the new state by the uni- 
tary operator, thereby constructing a new solution of the constraint. This proce- 
dure is followed for the two Dirac operators corresponding to fi((fi,Ptf>,p) '■= p<j> and 
f2{4>iP<t>iP) := P exp{=|=A/2K/3(</> — 0o)}- Notice that at <fi = 0o, these become the 
functions: p^ and p respectively and they act on the wavefunctions as the differential 
operator —ihd^ and as the multiplicative operator p, respectively. Explicit construction 
is as follows. 

Consider initial data, ^/±(fi, (p ), with the corresponding solution is denoted by, 

:= e^-^tf+^o) + e-^-<W\I/_Gu,0o) • 
Generate new initial data via the actions of p^, as, 

tUo#±CMo) := M^±(/i,0 o ) , p^±(fM,(po) ■= hV&V±(n,<l>o) . (4.35) 

Evolve these respectively, by e ±4 V / ^(<? i -</ , o)_ gy construction, these are solutions of the 
constraints, being of the form of (I4.33p . Explicitly, 

■= e^(^o)(nv^)vI/ + (^,0 o ) + e- i ^-^)(-^V / §)vI/_( /U ,0 o ) 
= -ihd^di,^) , (4.36) 

|/2U,flW) : = e^^H^+^^ + e-^^Vl*-^^) (4-37) 

Expectation values and uncertainties of these operators are used to track the quantum 
'evolution'. 

4. Physical inner product: 

It follows that the Dirac operators defined on the space of solutions are self-adjoint if 
we define a physical inner product on the space of solutions as: 

W) phys := " I duBijif (4.38) 

J 4>=4>o 

Thus the eigenvalues of the inverse volume operator crucially enter the definition of the 
physical inner product. For Schrodinger quantization, the integral is really an integral 
while for LQC it is actually a sum over \l taking values in a lattice. The inner product 
is independent of the choice of O - 

A complete set of physical operators and physical inner product has now been specified 
and physical questions can be phrased in terms of (physical) expectation values of 
functions of these operators. 

5. Semi- classical states: 

To discuss semi-classical regime, typically one defines semi-classical states: physical 
states such that a chosen set of self-adjoint operators have specified expectation values 
with uncertainties bounded by specified tolerances. A natural choice of operators for 
us are the two Dirac operators defined above. 

To be definite, let us consider the Wheeler-De Witt quantization. An arbitrary wave- 
function at some given </>o is expressed as an integral over k of the eigenfunctions Ck{p) 
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multiplied by some function of k and by a phase factor e iaj ( fc ^°. The inner product 
involves integral over /i. We can arrange to have a peak at some particular value 
by introducing a factor of — -v/3/2k|//*| so that for large k, the /x— integral will be dom- 
inated by /i ~ /i*. A large value of can be picked- up by choosing the function of k 
to be a suitable Gaussian. This particular large value can be related to a desired p*^. 
Here are the expressions: 

*semiGu,0o) := [dke-^^eMe 1 ^ -^ (4.39) 



k* = -vW^'Vi), , 0* = 0o- y/3/2K£n\fj.*\ . (4.40) 
The Gaussian allows the integrand to be approximated by 

efc „(^ e Mfc*)(>0-9!>.) ^ e ite|^|fc*+i(w(fc*) v /3/2K)to|/ i »,| ^ e O:*^n|/i/At*| 



where, in the last equality we used: k* is large and u{k*) ~ — a/2k/3 A;* (see the 
footnote [TO]) . The integral in the inner product will now pick-up contribution from 
near fi ~ /i*. With these observations, it is easy to verify that the 'initial' semiclassical 
wave function given above gives (p^) = p*i and (|/i|</, ) = /V The initial semiclassical 
wavefunction evolves into ^/ se mi(A i , 0) which is same as the initial wave function with 

00^00 

For LQC, the e^p) functions are different [T7] and the physical expectation values are 
to be evaluated using the physical inner product defined in the LQC context. 

6. Evolution of physical quantities: 

We have now the physical wave function, evolved from \I/ S emm- Since it retains the form 
of the initial wavefunction, the k integral can be approximated as before and thus will 
lead to same expectation value for p^ for all <fi. For the expectation value of (|/^|^ o )(0), 
the fi integral will be saturated by the new phase, — O + a/3/2k £n\ /////* | — 0. And 
this we recognize as precisely the solution ( 14. 7p . 

Thus, in the WdW quantization, the classical relational evolution, p = p((f>) is repro- 
duced by the expectation values of the Dirac observables. 

Exercise: compute/estimate the uncertainties of the Dirac observables in the semiclas- 
sical state given above.. 

7. Resolution of Big Bang Singularity: 

A classical solution is obtained as a curve in (//, 4>) plane, different curves being labelled 
by the points (/i*,0*) in the plane. The curves are independent of the constant value 
of These curves are already given in (14 .7p . 

Quantum mechanically, we first select a semi-classical solution, ^/ semi (p^, \x* : </>) in 
which the expectation values of the Dirac operators, at (p = (fi , are p^ and fi* respec- 
tively. These values serve as labels for the semi-classical solution. The former one 
continues to be p^ for all (ft whereas (|a*|^ o }(0) =: Mp*,^* (0), determines a curve in the 



n In the above heuristic reasoning for the form of the semiclassical state, we have glossed over some 
technical issues such as whether the states exhibited are in the domain of the O operator which requires 
carefully stipulating conditions on the function of k (which has been taken to be a Gaussian). These are too 
technical to go into here and are not expected to affect the corresponding discussion for LQC. A discussion 
of these issues may be seen in [15] . 



43 



(p, (j)) plane. (We determined this curve above, using stationarity of the phase for the 
Schrodinger quantization). In general one expects this curve to be different from the 
classical curve in the region of small p (small volume). This is what happens for the 
loop quantized theory. 

The result of the computations is that Schrodinger quantization, the curve |^/| p * }/x * (</>), 
does approach the fx = axis asymptotically. However for LQC, the curve bounces away 
from the fx = axis. In this sense - and now inferred in terms of physical quantities 
- the Big Bang singularity is resolved in LQC. It also turns out that for large enough 
values of p^, the quantum trajectories constructed by the above procedure are well 
approximated by the trajectories by the effective Hamiltonian. All these statements 
are for semi-classical solutions which are peaked at large p* at late times. 

Two further features are noteworthy as they corroborate the suggestions from the effective 
Hamiltonian analysis. 

First one is revealed by computing expectation value of the matter density operator, p ma tter := 
§fc) 2 |p|~ 3 , at the bounce value of \p\. It turns out that this value is sensitive to the value 
of p*^ and can be made arbitrarily small by choosing p*^ to be large. Physically this is 
unsatisfactory as quantum effects are not expected to be significant for matter density very 
small compared to the Planck density. This is traced to the quantization of the gravitational 
Hamiltonian, in particular to the step which introduces the ambiguity parameter p . A novel 
solution proposed in the "improved quantization", removes this undesirable feature. 

The second one refers to the role of quantum modifications in the gravitational Hamiltonian 
compared to those in the matter Hamiltonian (the inverse volume modification or B(p)). 
The former is much more significant than the latter. So much so, that even if one uses the 
B(fi) from the Schrodinger quantization (i.e. switch-off the inverse volume modifications), 
one still obtains the bounce. So bounce is seen as the consequence of 9 being different and 
as far as qualitative singularity resolution is concerned, the inverse volume modifications are 
un-important. As the effective picture (for symmetric constraint) showed, the bounce occurs 
in the classical region (for j = 1/2) where the inverse volume corrections can be neglected. 
For an exact model which seeks to understand why the bounces are seen, please see [T9] . 

Improved Quantization: The undesirable features of the bounce coming from the classical 
region, can be seen readily using the effective Hamiltonian, as remarked earlier. To see the 
effects of modifications from the gravitational Hamiltonian, choose j = 1/2 and consider the 
Friedmann equation derived from the effective Hamiltonian leading to the effective energy 
density ( 14.26[) . with matter Hamiltonian given by H maittev = ^p^\p\^ 3 ^ 2 - The positivity of 
the effective density implies that p > p* with p* determined by vanishing of the effective 

energy density: p* := p c i(p*) = ( 87rG ^" 7 p*)' 1 . This leads to = ^/ 47rG H 7 |^| and 

p* = \/2( 87rG g ° 7 )~ 3 / 2 \p ( f > \~ 1 . One sees that for large |p^|, the bounce scale |p*| can be large 
and the maximum density - density at bounce - could be small. Thus, within the model, 
there exist a possibility of seeing quantum effects (bounce) even when neither the energy 
density nor the bounce scale are comparable to the corresponding Planck quantities and this 
is an undesirable feature of the model. This feature is independent of factor ordering as long 
as the bounce occurs in the classical regime. 

One may notice that if we replace po — > p(p) := y/A/\p\ where A is a constant, then the 
effective density vanishes when p c i equals the critical value p CI1 t '■— 
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independent of matter Hamiltonian. The bounce scale is determined by p* = p cr i t which 
gives | p ,| = (_£L)V3. Now although the bounce scale can again be large depellding uponp0 , 

■^Pcrit ^ 

the density at bounce is always the universal value determined by A. This is a rather nice 
feature in that quantum geometry effects are revealed when matter density (which couples 
to gravity) reaches a universal, critical value regardless of the dynamical variables describing 
matter. For a suitable choice of A one can ensure that a bounce always happens when the 
energy density becomes comparable to the Planck density. In this manner, one can retain 
the good feature (bounce) even for j — 1/2 thus "effectively fixing" an ambiguity parameter 
and also trade another ambiguity parameter p for A. This is precisely what is achieved by 
the "improved quantization" of the gravitational Hamiltonian |20j . 



The place where the quantization procedure is modified is when one expresses the curvature 
in terms of the holonomies along a loop around a "plaquette". One shrinks the plaquette 
in the limiting procedure. One now makes an important departure: the plaquette should 
be shrunk only till the physical area (as distinct from a fiducial one) reaches its minimum 
possible value which is given by the area gap in the known spectrum of area operator in 
quantum geometry: A = 2y/3TcyGh. Since the plaquette is a square of fiducial length /i , its 
physical area is pl\p\ and this should set be to A. Since \p\ is a dynamical variable, po cannot 
be a constant and is to be thought of a function on the phase space, p(p) := i/A/|p|. Thus 
we need to define an operator corresponding to the classical expression: hf := exp(i|/(p)c), 
we have taken a general function f(p). This is little non-trivial since there is no c operator 
and c, p are conjugate variables. 

Observe that the usual holonomy operator effects a shift in the argument of eigenstates of 
the triad operator and formally the operator looks like exp(u-^) i.e. it effects the action of 
a finite diffeomorphism generated by a vector field on the wavefunction. We will take this 
as a guiding principle. 

Let denote a diffeomorphism effecting a unit parameter shift along the integral curve 
of the vector field /(/■*) and the corresponding pull-back map. We define hf^(p) := 
[$*(\]/)] : = \I>($j(p)). As argued above, for a constant function, this reduces to the usual 
action (14. 9p . It can be checked directly that this action is also unitary in the kinematical 
Hilbert space: ($, \&) = ^2 ^*(t i )^(f i ) where the sum is over a countable set (this follows 

from := £ Mecountable subsetcR ^(/i)|p)0 



To compute a unit parameter shift due to the diffeomorphism generated by f(p), solve the 
equation 

rv+l 

dv = 1 (4.41) 



dx 



/(*) 



This will give p! := <&f(p). 

For the specific choice of f(p) := p(p) := v^A|p| -1 ^ 2 , A := 7a/3£|/4, p(p) = *y£pp/6, one 
gets, 

/ A I O . /O 





sgn(//)l//| 3/2 = sgn(p)\p\^ + K- 1 K]= 3V^ (442) 




2 For a comparison in Schrodinger quantization, see remarks in |17j . 



45 



(A*) 



It is evident from the above that if we define v := i^sgn(/i)|/i| 3//2 , then the middle eqn reads: 
v' — v + 1. This suggests that we use \v) as a basis instead of Apart from the constant 
K, and the sgn, v is related to the eigenvalue of the volume operator Ipl 3 / 2 . Note that v 
as a function of fi is one-to-one and on-to. Using the hp operator and a basis labelled by 
volume eigenvalues, the Hamiltonian constraint is defined and difference equation is obtained 
as before. The relevant expressions are: 



V\v) 



\p\- 1/2 



3=1/2,1=3/4 

B(v) 

©Improved^ (v,<j)) 



iv) 



Ksgn(n)\n\ 



3/2 



7\ 3 / 2 4 I I I \ 



(4.43) 
(4.44) 



*(u + k) 

3 (j§; 

2 V 6 

3N 3/2 

2 



-1/2 



(4.45) 

K 1/3 \v\ 1/3 \\v + 1| 1/3 -\v- 1| 1/3 | *(v)(4.46) 

1/3 | 3 



K\v\ \\v + 1| 1/3 - \v 

3ttKG 



(4.47) 



>) + C°(v)V(v, 
C-(v)ty(v-4, 



\v + 2\ | |v + 1| - |« + 3| 
-4) , C°(^) := -C + 



(4.48) 
(4.49) 
(4.50) 



Thus the main changes in the quantization of the Hamiltonian constraint are: (1) replace 
pLo — ?• fb := a/A/|p| in the holonomies; (2) choose symmetric ordering for the gravitational 
constraint; and (3) choose j = 1/2 in both gravitational Hamiltonian and the matter Hamil- 
tonian (in the definition of inverse powers of triad operator). The "improvement" refers 
to the first point. This model is singularity free at the level of the fundamental constraint 
equation (even though the leading coefficients of the difference equation do vanish, because 
the parity symmetry again saves the day); the densities continue to be bounded above - 
and now with a bound independent of matter parameters; the effective picture continues to 
be singularity free and with undesirable features removed and the classical Big Bang being 
replaced by a quantum bounce is established in terms of physical quantities. 

There is yet another spin on the story of singularity resolution! 



4.2.1 Madhavan Quantization [18]: 

The improved quantization scheme works primarily through the holonomy corrections, so 
much so that even if the inverse volume corrections in the matter are turned-off by hand, 
the singularity resolution continues. Madhavan works within the same kinematical Hilbert 
space of LQC but treats the Hamiltonian constraint differently, exploiting its specific, simple 
form for the massless scalar matter. In his quantization of the Hamiltonian constraint, it is 
the inverse volume corrections that are responsible for singularity resolution (also in terms 
of physical quantities) and holonomy corrections are by-passed completely. 
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He observes that the classical Hamiltonian constraint (I4.27p . is quadratic in c and is of the 
form of difference of two squares. It can therefore be written as a product, 



tot — —C + C- 



C± :=- 




The C± are linear in c. Since to define physical Hilbert space, a general procedure is to 
average over the group generated by the constraint, and this involves exponentiation of the 
constraint, one can directly define these operators which involve factors of the same form as 
the hf operators of the improved quantization! 

The key differences in Madhavan quantization are: (1) The Hamiltonian constraint is reg- 
ulated differently from the analogue of LQG and as a consequence, there are no holonomy 
corrections (sin 2 (c)). However, the inverse triad corrections can be incorporated in the defi- 
nition of C± through the p^ term. One again uses the volume eigenvalues basis for the inverse 
triad definition; (2) the physical states are constructed directly by group averaging using the 
well-defined unitary operators, e iaC± , consequently there is no difference/differential equation 
to be solved; (3) One of the Dirac observables, p^ is the same but another one is somewhat 
different. Nevertheless, classical solutions can be derived from their expectation values; (4) 
The issue of independence from the fiducial cell (discussed in the next subsection) is also 
addressed differently. 

The results are: (i) without inverse triad modifications, the classical (singular) solution is 
recovered; (ii) with inverse triad modifications, there is extension of the solution past the 
classical singularity with the energy density remaining bounded all through, making the 
extension non-singular. There is no bounce, but a regular extension! 

Although Madhavan's procedure of bypassing the holonomy corrections completely, is tied 
to the particular form of the constraint of the isotropic model (and hence may not ex- 
tend to other models), it does demonstrate the possibility that there are inequivalent ways 
of constructing physical Hilbert space and observables starting from the same kinematical 
structures. Secondly, singularity resolution need not be seen only as a classical/quantum 
bounce, a regular extension is also a distinct possibility. 

More details should be seen in [18]. 

4.2.2 Role of the Fiducial cell in spatially flat models 

Recall that in the description of spatially homogeneous and isotropic models one begins with 
a metric of the form ( |4.ip . The spatial metric is a metric with spatially constant (but possibly 
time dependent) curvature. This is conveniently taken to be a time dependent scaling of 
a fixed co-moving metric with corresponding co-moving coordinates. Although not strictly 
necessary, let us assign length dimension to the co-moving coordinates and take the scale 
factor to be dimensionless. For non-flat models the co-moving metric can be normalized to 
have the Ricci scalar to be ±1 in appropriate units (Ricci scalar has dimensions of (length) -2 ). 
Note that this is a local condition, and by homogeneity, holds everywhere on the spatial 
manifold. It is independent of the size of the spatial manifold. For flat models, such a 
normalization of the co-moving metric is not possible. In this case, there is an arbitrariness 
in the definition of the scale factor. Clearly, by focusing only on those quantities which are 
invariant under constant scaling of the scale factor, eg a/a,d/a, the energy density etc we 
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can obviate the need for choosing a co- moving metric/coordinates. The equations of motion 
- the Friedmann equation and the Ray chaudhuri /continuity equation - reflect this feature. 

However, spatial flatness, homogeneity and isotropy also implies existence of (global) Carte- 
sian coordinates with a metric g^ = 8ij with the coordinate differences giving distances in the 
chosen unit of length. This unit is arbitrary, but also determines the unit of time by putting 
speed of light to be one. Change in this unit results in an overall scaling of the space-time 
metric but does not affect the scale factor. The scale factor is now unambiguously identified 
and co-moving coordinates and metric are also fixed. 

Construction of a quantum theory of the scale factor degree of freedom (and matter homo- 
geneous degrees of freedom) begins with a four dimensional action principle restricted to 
homogeneous modes of the fields. The action contains a spatial integration which is diver- 
gent for spatially flat models, thanks to homogeneity. To have a well defined phase space 
formulation, we need to regulate this divergence. This is done by introducing an arbitrarily 
chosen fiducial cell, specified by finite ranges of the co-moving coordinates (thus having a 
finite co-moving volume Vq) and restricting the integrations to this cell. Note that this need 
and the freedom in the choice of the fiducial cell arises strictly due to the need for an action 
formulation for the full theory and the assumption of spatial homo geneityf^l. All subsequent 
computations will carry a dependence on this cell, either explicitly or implicitly. In the end, 
this dependence is to be removed by taking a suitable limit V — > Precisely at what 

stage and how should one take the limit? 

In the canonical formulation of the full theory, the fiducial volume, Vo, appears in the sym- 
plectic structure. This can however be absorbed away by redefining the canonical coordinates 
(c,p — > c,p). This makes the canonical coordinate p to have dimensions of (length) 2 . Note 
that the physical volume of the cell is a 3 (t) x Vo is now directly given by |p| 3 ^ 2 . Apparently, 
there is no reference to the fiducial cell any more in the model. However this is not so. 
The |p| 3//2 is the physical volume of the fiducial cell. All subsequent computations, whether 
classical or quantum, done using the (c, p) variables^! have no explicit reference to the cell. 
For example, the classical solution obtained in terms of phase space trajectory, (eqn. I4.7p . 
does not depend on Vo. 

As discussed above, the subsequent steps in the quantization, do not introduce any further 
dependence on the fiducial cell. It is no where in sight even in the computation of the phase 
space trajectory (expectation values of the Dirac observables) . These trajectories of course 
differ from the corresponding classical trajectories. The problem of Big Bang singularity is 
however phrased in the framework of space-time geometry, specifically, in terms of backward 
evolution of the scale factor. So we need to transcribe the phase space trajectories (computed 
in terms of expectation values of Dirac observables) into evolution of the space-time geometry 
i.e. the scale factor. At this stage, a scale factor (and an explicit reference to the fiducial 
cell) is re-introduced via the triad variable as, a := ^ A /p where £ has dimensions of (length)" 1 

and can be identified with the fiducial volume: = Vq^ 3 (since p 3//2 is the physical volume 
of the fiducial cell). The phase space evolution then gives a(t). The scale factor evolution 

13 The action formulation (Lagrangian or Hamiltonian) in turn is required for a quantum theory. The 
classical theory needs only equations of motion which are independent of any cell. 

1 The limit Vo — >■ oo can be viewed as a convenient way to pick-out Vo— independent terms and/or could 
also be heuristically motivated by noting that the definition of homogeneity identifies the spatial manifold 
with the group manifold and this group manifold is R 3 for the present case. 

15 During the process of loop quantization, fiducial scales could appear again eg through the holonomies 
along edges. However as explained in the footnote EJ the Vq disappear. 
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so deduced could have some dependence on £. After taking the limit £ — > (Vo — > oo), 
the evolution that survives, is the prediction of the quantum theory, Whether or not this 
evolution is singularity free (i.e. all physical quantities remain bounded through out the 
evolution) is the central question of interest. Since the phase space curves are inferred from 
the expectation values, the states in which these are computed are also important to specify. 
These are expected to be computed in physical states peaked on large volume p>fp and 
small energy density (say), corresponding to a classical regime. The singularity free evolution 
is required to hold for all such states. 

The classical evolution given in eqn (14. 7p provides an example of this transcription. We see 
that £ -2 cancels out from both p and p* and the classical evolution of the scale factor is 
independent of the fiducial cell as it should be. The LQC computed solution for p, always 
shows a bounce, is not very explicitly expressed and also contains an implicit dependence on 
the fiducial cell. It matches pretty closely with the classical solution in the large p regime 
and therefore could be expected to be Vo independent in these regimes. This removes the 
cell dependence in the initial condition and the question boils down to whether the bounce 
feature and value of p at the bounce, is independent of V . 

The APS investigations [TU| [TT] found that, in the /j,q— scheme, a bounce can occur even for 
low values of energy density something which is not exhibited by the observed isotropic 
universe. Furthermore, the energy density at the bounce - which is a physical observable - 
has a Vo— dependence. So the quantization scheme has some problems. What exactly does 
this mean? 

Note that this does not necessarily mean that there is any mathematical inconsistency in 
the process of quantization. However, the constructed quantum theory should agree with 
GR for low energy densities (i.e. have acceptable infra-red behaviour) and hopefully also 
imply a non-singular evolution. It is possible that it may fail this expectation. APS analysis 
concludes that the fiQ— scheme with symmetric ordering fails this test. 

In retrospect, this failure could have been inferred in the following way. The earlier methods 
of analysis were based on WKB approximation and effective Hamiltonians|12j derived from 
it. This allowed us to encode quantum modifications in terms of effective density and effective 
pressure, defined by computing the left hand sides of the Friedmann and the Raychaudhuri 
equations using the effective Hamiltonian. In these papers, the scale factor was introduced 
by setting p := a 2 /4 (dimensionful) and therefore still refers (implicitly) to fiducial cell. To 
make the fiducial cell explicit, replace this scale factor as a — > a^" 1 . The expressions for 
the energy density (say) can be transcribed in terms of the (dimensionless) scale factor and 
£. As explained above, prediction of the quantum theory for the scale factor evolution is 
obtained by taking the limit V — > oo(£ — > 0). Note that this is now done at the level of 
equations as opposed to at the level of individual curves which need initial conditions to be 
chosen (which we have argued to be independent of Vo). 

Recall the effective density given in eqn. |4.26] . 



matter 




where 



(4.51) 
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The square bracket contains the modification implied by holonomy correctional The inverse 
volume corrections are contained in the matter Hamiltonian (see eqns. (14.121 14. 13j) ) . 

Consider first the inverse volume corrections. Observe that with p — > a 2 £~ 2 /4 , q ~ a 2 £~ 2 ^> 
1 in the limit £ — > 0. The limiting form of Fg then implies that p~ 3 / 2 -£f matter ~ v\^ a ~^ ■ F° r 
massless scalar matter, the classical equation of state has P/p = 1 and hence the classical 
density behaves as ~ a -6 . Thus p 2 £ 6 must be a constant for any particular solution. For 
the leading term to give the classical evolution, we have to take the limit £ — > along with 
p^ — > oo keeping p^ 3 a constant specifying a particular initial condition. This understood, 
the p~ 3 / 2 -£^matter factors go over to the cell independent classical density plus corrections down 
by q~ 2 in the limit. Thus, inverse volume corrections, simply vanish when Vo — > oo is 
imposed. 

Now consider the holonomy corrections. By the same logic as above, the holonomy correc- 
tions, second term in the square bracket in the effective density expression, goes as £ -2 — > oo!. 
This is clearly unacceptable. Thus, in the Hq scheme of quantization, the inverse volume mod- 
ifications do not survive the limit while the holonomy modifications give an inconsistency 
and neither shed any light on the singularity resolution issue. 

True as these features are, they are not immediately conclusive to look for alternative quan- 
tizations because the fault may be with the WKB approximation and the corresponding 
effective Hamiltonians. For instance if one took the effective density from the first paper 
of [12], the holonomy corrections would also be down by inverse powers of q and would 
vanish which is okay for the classical limit but the extrapolation to the quantum regime is 
unreliable since WKB is unreliable at turning points. Perhaps, physical level computations 
would clarify the issue. This is indeed the case. With the physical level computations, APS 
results show the unacceptability of the /i — scheme while in Madhavan's approach, with no 
holonomy corrections, the inverse volume corrections would simply vanish by the argument 
given above. Both APS and Madhavan have also suggested ways out. 

The APS analysis discussed above shows that the /x — > fi(p) = V~K£p/y/p substitution in 
the holonomies used in replacing the c variable, suffices to obtain a non-singular evolution 
with good infra-red behaviour. It implies that the deviations from classical evolution (eg 
close to the bounce) occur when the energy density reaches a universal, maximal value. 
This substitution also renders the inverse triad correction from the matter sector highly 
suppressed^!. 

Madhavan suggests that along with the APS suggested substitution, one should also intro- 
duce a multiplicative parameter A as /i — > A/i(p). Since only corrections that drive the quan- 
tum modifications are the inverse volume correction and these go as (X/vi) 2 ~ (A£ 3 £p/a 3 ) 2 
(see eqn 15.451) . £ independence is achieved by choosing A ~ (£/£) 3 where the new dimen- 
sionful parameter £ is supposed to reflect a scale provided by an underlying LQG state 
supporting the homogeneity approximation. For details, please see [18J. 

To summarize: For the spatially flat, isotropic models, let us choose the Cartesian coordinates 
with the standard Euclidean metric for the spatial slice and choose proper time as the time 
coordinate so that the space-time metric takes the form (14. ip . The scale factor is now 
specified unambiguously. For constructing a quantum theory, we need to choose a regulator 

16 In the first paper of [12], the effective density contained only the leading terms of the holonomy corrections 
which have been summed up in the second paper. 

17 From eqn. 15.451 with A = 1, one sees that the corrections go as v± ~ £pp~ 3 ~ £p(£) 6 a~ 6 — > as £ — >• 0. 
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cell with volume Vq. While the cell dependence can be hidden by choosing scaled variables, 
it manifests again because quantum computed evolution must be transcribed in terms of 
the scale factor evolution. This is necessary because the classical Big Bang singularity is 
understood as a singular evolution of the scale factor so its resolution lies in making the 
evolution non- singular. 

The scale factor evolution can be cast in the form of the Friedmann equation with possible 
deviations from classical evolution, encoded in the effective density. A prediction of the 
quantum theory is the surviving correction terms after taking the limit Vq — > oo. A quan- 
tum theory could be understood to have resolved the Big Bang singularity if the surviving 
evolution is non-singular. Cell independence of quantum corrections automatically implies 
that non-trivial limit exists. Not every quantization scheme passes this test. 

There are two types of corrections - the holonomy corrections and the inverse triad correc- 
tions. These have different properties in the limit. The APS quantization with lxq— scheme 
implies that holonomy corrections dominate and lead to unphysical implications. These are 
cured by the Jx— scheme. The Madhavan quantization scheme, even with the p, substitution 
in the inverse volume definition, these corrections again vanish unless additional A param- 
eter is introduced. In either case, extra ingredients (scales) have to be 'imported' to get 
non-trivial results. Both the schemes have ingredients (role played by the area operator in 
the APS schemed and the specific form of constraint in the Madhavan scheme) which do 
not have counterparts in the full theory as it is understood at present. 

So far the discussion has been within the context of full theory being classically reduced 
directly to a homogeneous and isotropic model. In the next subsection, we briefly discuss how 
homogeneity and isotropy can be viewed from within a particular quantized inhomogeneous 
model. 



4.2.3 A View from Inhomogeneity 

To keep the flow of the in focus, basic details of the inhomogeneous model are given in the 
appendix 5.4. 

The fundamental change in the way homogeneity is viewed, is that it is a property exhibited 
by a state of an inherently (spatially) inhomogeneous model. For definiteness, a lattice model 
with a lattice spacing £ is taken. This allows for states of the model which can be considered 
as homogeneous on a certain scale, eg IqN 1 ^ 3 . This also allows the fields to be restricted to 
be periodic on this scale. Thus what is fundamental is the lattice spacing Iq below which it 
makes no sense to consider inhomogeneity (or inhomogeneities are not probed) and a scale 
A rl / 3 ^o provided by a state of the model. Let us call the former as the micro-scale and the 
latter as the macro-scale. The fundamentally isotropic model refers to the macro-scale. The 
fiducial cell of the isotropic model is determined by these two scales with fiducial volume 
given by, Vo := N£q. Notice that in this view, Vq (or N) is a property of a quantum state 
and there is no reason to contemplate a limit Vq — > oo. 

Now all quantum effects due to inverse volume and holonomies, arise at the micro-scale. 
To see how these translate or correspond to the quantum effects seen in the fundamentally 

18 The logic used to motivate the role of area operator, is also extended to other Bianchi models such that 
when isotropy is imposed, the /2-scheme of isotropic model is recovered backpTj. 



51 



isotropic model, let us begin by identifying variables. 

In the lattice model, isotropic connection is defined by k(x) = c, V x, I. Identifying this 
constant value with Ci SO and comparing the basic link holonomies with the holonomies of the 
isotropic model implies, 

c iso := V 1/3 c iso = V 1/3 ~kj := V^tfc^ = iV 1/3 c lat (4.52) 

The first equality is the definition from the isotropic model, the second one identifies the 
isotropic connection with the lattice connection k, the third one defines c\ at and forth equality 
gives the final relation between the 'c' variables of the isotropic and the lattice models. 

The state exhibiting isotropy on the macro-scale, may be characterised by stipulating that 
the pi(v) values are all mutually equal, and equal to pi at , for the vertices comprising the 
fiducial cell and that this value is identified with the isotropic variable, p; so . Recall that the 
isotropic operator is obtained by an averaging of the lattice operators. The stipulation says 
that the average value is realized at each of the lattice vertices. This leads to, 

p iso := V* /3 p iso = V 2/3 put = V 2/3 put£o 2 = PiatiV 2/3 (4.53) 

These identification give a relation between isotropic variables and a state of lattice model 
with a scale parameter N. A dynamically evolving isotropic universe may be thought of as a 
family of lattice states with the scale N being a function of the volume eg larger number of 
elementary cells get 'homogeneised/isotropised'. As an example, if N oc p 3 ^ 2 , then pi at will 
be a constant! The expressions seen in the context of isotropic models with fiducial cell of 
size Vo are now to be applied to the elementary cell of size £q. 

In view of these identifications, let us consider the two specific corrections seen in the isotropic 
model, namely the the inverse volume corrections and the holonomy corrections. In the 
lattice model, these corrections arise in the same manner as in the isotropic model, but 
from the micro-cell. The inverse volume corrections are in powers of p*/p = p*/pi&t = 
{p*/Piso)N 2 / 3 . For N oc pf^ 2 , these are independent of pi SQ . (ii) With the holonomy corrections 
included, the effective density (I4.5ip . is of the form p(l — p/pcnt) with p~ it ~ ^7 2 pi a t = 
k^ 2 P\s N~ 2 / z . Again for iV cx pf^ 2 , p cr i t is independent of pi SD . For a different dependence of 
iV on p iso , the corrections will have non-trivial dependence on p iso and therefore on Vo, but Vo 
is no longer a purely mathematical artifact but is dictated by the underlying inhomogeneous 
state. 

In effect, a perspective from an underlying inhomogeneous model suggests that the fiducial 
cell of a homogeneous model is selected by a state of the inhomogeneous model and a 
dependence on the Vo := £qN is not necessarily unphysical. For a more detailed discussion 
of ramifications of these ideas, please see [16J. 
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Chapter 5 
Appendix 



5.1 Symmetric connections 

We assume that we have a manifold M on which are defined connection A^(z) with a taking 
values in the Lie algebra G of a gauge group G. Assume further that there is an action of 
a symmetry group S on M under which we want to have appropriate notion of invariance. 
The infinitesimal action of the symmetry group is generated by a set of vector fields 
which represent the Lie algebra of S: [£ m , £ n ] = f p mn X p . The Lie algebra of G is generated 
by matrices T a satisfying: [T b ,T c ] = C a bc T c . 

When we have an ordinary tensor field, T, on a manifold, it is defined to be invariant under 
(or symmetric w.r.t.) the action of an infinitesimal diffeomorphism generated by a vector 
field £ if its Lie derivative with respect to £ vanishes: L(T = 0. When the tensor fields also 
transform under the action of a gauge group, then the invariance condition allows the Lie 
derivative to be an infinitesimal gauge transform of the tensor field: L(T = 6w(£)T, where 
W(£) is valued in the Lie algebra G. Notice that this associates a gauge transformation with 
a diffeomorphism. 

This association has to satisfy two conditions: (a) If we took a gauge transform of the tensor 
field and then applied the diffeomorphism, the defining condition must be gauge covariant: 
L^(T 9 ) = 8w9(£)T, where W 9 (£) = g~ 1 W(£)g + g~ 1 L^(g), and (b) The Lie derivatives repre- 
sent the Lie algebra of the symmetry group: [L u , L in \T = L [u ^ n] T = f p mn L^ p T and the W m 
must obey the consequent conditions. The task is to find those tensor fields which satisfy 
the invariance conditions subject to the allowed gauge transformations. 

This is aided by another consequence of the symmetry action. The action of the the symmetry 
group S on M implies that M can be expressed as a collection of orbits of S. We will assume 
the simpler case where all orbits are mutually diffeomorphic and are given by S/ F where F 
is the stability subgroup of the S— action. Thus we obtain M ~ B + S/F. Here S/F is an 
orbit which is necessarily a coset space while B is a manifold whose points label the orbits. 
Note that a non-trivial subgroup F of S means that a subset of vector fields £ m vanish at 
some point. Corresponding to this structure of M, its tangent and cotangent spaces are also 
decomposed. 

The solutions of the invariance conditions are constructed by using the available structure 
on the group manifold S and projecting these onto S/F. Here are some details for the gauge 
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connection^ [22] . 

Let A := A a ^T a dx^ denote the G valued connection 1-form (the gauge potential). Un- 
der a gauge transformation it transforms as: A 9 := g~ 1 Ag + g~ l dg and infinitesimally, 
g = 1 + eW, 5 eW A = eD(W) : = e(dW + [A,W\), W := W a T a . Under an infinitesimal diffeo- 
morphism, x M := x M + e£ M , it transforms as 5 e ^A := —eL^A = —t{d^ u A u + ^ u d u A^)dx^ = 
- {i^dA + d(i^A)}. 

The invariance conditions, when there are many symmetries, are: 

L im A = D(A) W m := dW m + [A, W m ] & d£ v A v + = d,W m + [A,,, W m ] . (5.1) 

Here, W m are some G valued scalars on the manifold M associated with £ m . We induce a 
gauge transformation on the W m by demanding that the above condition be gauge covariant: 

L U A 9 := D(A°) W* ^ := g- l W m g + g-'L^g , L u g := ^ d, g (5.2) 

The Lie algebra of the vector fields, £ m , implies that the W m 's must satisfy: 

[L im ,L in ]A = L [UM A = f p mn L^ p A => 

D(W mn ) = F rnn D(W p ) , W mn :=L im W n -L in W m + [W m ,W n ] , or 

= D(PU mn - ^ mn PU p ) (5.3) 

where [W m , W n ] is the bracket in the Lie algebra, G. 

Exercise: For the field strength F^ u T a , verify that L^ m F = [F,W m ]. 

Exercise: Let E be a vector field valued in G which transforms as E 9 = g~ 1 Eg. The condition 
for symmetric E would be L^ m E = [E,W m ). Show that the W m transforms as before and 
the symmetry Lie algebra implies [E, W mn — f p mn W p ] = 0. 

Suppose that Xmn '■= W mn — f p mn W p ^ 0. Then the above equations imply conditions on 
symmetric field strength (and indeed on all symmetric quantities). For example, Dx mn — 
implies that [F, Xmn] = Qj- This would mean that the field strengths must commute with 
the Xmn's. We will assume that there does not exist any x valued in G such that Dx = 0. 
This implies that Xmn = i.e. 

L u W n - L u W m + [W m , W n ) - f p mn W p = . (5.4) 

Note that this is a condition involving only the W„'s and the symmetry generators £ m 's. Also 
observe that for vector fields corresponding to the stability subgroup F, the (15. 4 p reduces to 
[W m , W n ] = f p mn W p at the points where these vector fields vanish. 

The task is to characterise the symmetric connections satisfying eqn. ( 15. ip with W n satisfying 
eqn. fl5.4[) modulo gauge transformations ( 15. 2 p on a manifold M ~ B x S/F. The strategy 
is to show that the gauge freedom allows W n to be taken in appropriate form and then 
determine the form of symmetric connections in the same gauge. 

the context of Kaluza-Klein approach to unification, the analysis of invariant quantities was carried 
out to construct suitable ansatz. There, the space-time is taken to be of the product form M4 x B with B 
a compact manifold. The isometry groups of the compact manifold played the role of symmetries and the 
forms of the fields on the space-times were obtained. 

2 Note that D\ = follows only for the symmetric connections. If it were to hold for all connections, then 
[F, x] = would hold for all field strengths and this would correspond to a reduction of the gauge group to 
the little group of x, nn . This is analogous to the non-trivial Higgs vacua situation. 



54 



Consider first the case where F = {e} so that S/F = S itself. Introduce local coordinates 
(x l ,y a ) on B x S. Without loss of generality, we can take the symmetry generators to 
be functions of y and with zero components along B. Noting that the vector fields £ m on 
S are independent, it follows that the matrices ^(y) are invertible and therefore we can 
define new G valued 1-forms as: W m (x,y) := ^W a (x,y). It is easy to see that (i) W a 
transform exactly as a G— connection and (ii) the condition ( 15. 4p is just the statement that 
this connection is flat. Therefore, locally it is always possible to choose W a (x,y) = and 
hence W m (x, y) = 0. Having chosen W m 's to be zero, the gauge transformation freedom 
is restricted to L^ m (g) = ^d a g = i.e. the gauge functions must depend only on the x 
coordinates. 

In this gauge, the invariance conditions can be written separately for /j, = i and for // = a 
as: 

^d a Ai = and (L u A) a = 

The first implies that A{ depends only on x and so do the gauge transformations. Hence Ai 
is a G— connection on B. The solution for A a is obtained as follows. 

On a group manifold, there are left and right actions of the group onto itself which commute. 
Consequently, these generate left (right) invariant vector fields and 1-forms (the Maurer- 
Cartan forms). Apply these to the group manifold S. Assume that £ m generate left action 
on the group manifold and use the left invariant 1-forms eg the unique, 5— valued Maurer- 
Cartan form 0mcB It is immediate that L^Qmc — 0. To obtain a G— connection, we need 
a map A : S_ — > G. Given such a map, we can define A := A(Q M c) A a a := ^(0^)™. 
Now L^ m A = = L^ m ($)©7v/c* + implies that the "Higgs" fields, $ m are constants on S 
i.e. are functions only of x l . 

Thus, for the case where the S— action is free (F is trivial), the symmetric connections can 
be written as A = Ai(x)dx % + <& m (x)uj m a {y)dy a where A is an arbitrary G— connection on B 
and $ m (x) are "Higgs" scalars valued in G. 

When the S— action is not free, the vector fields £^ are tangent to S/F and the above steps 
do not go through immediately. Nevertheless we can still find invariant, S_— valued 1-forms 
on S/F and a suitable map A to construct invariant connections. To see this, note that 
there is the natural projection map it : S —> S/F. Choose an embedding i : S/F — >• S. As 
discussed above, on S we have vector fields £ m generating left action and the corresponding 
Maurer-Cartan form Qmc- Using 7T* we push-forward the vector field on to S/F and using i* 
we pull-back Qmc on to S/F. The projected vector fields match with the £ m (by definition 
of the symmetry action). Thus we get, 

£n:=7r*(£n) , u:=i*{Q MC ) ; L^ n Q MC = L^u = . 

As before, u is valued in S_ since Qmc is- Introduce <3>^ as before and define A a a := 
Using these definitions let us rewrite the defining equations as: 

Invariance condition : ( L U®k ~ [$k, W m ]) = L^ n W m ; (5.5) 

Lie Algebra condition : L u W n - L in W m - [W m , W n ] = f p mn W p ; (5.6) 
Gauge transformations : = g~ x W m g + g~ x L^ m g ; (5.7) 

: {i> k a )H-9~^kg) = g-'L^g (5.8) 

3 For classical matrix groups, these are given by g~ 1 dg and dgg^ 1 and are left and right invariant respec- 
tively. 
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In writing the first equation we have used L^ m oj n = and also multiplied by 

The last equation implies that $ n 's transform as the adjoint representation of G iff the 
gauge transformations are constant over S/F. We can take $ n to transform by the adjoint 
representation of G, thereby restricting the gauge transformation to be constant over S/F. 
In such a case, W n 's also transform the same way. For trivial F, we can transform away W n 
to zero and recover the previous case. For non-trivial F, this is not the case. 

For non-trivial F, there exist a point, y say, in S/F at which the vector fields £ m ,m = 
1, • • • , dim(F) vanish. Then Lg m terms drop out. Consider the equation ( 15. 6ft for m,n. 
Then, at y , we must have \W m , = f~ m nW p . Since F is a subgroup, the sum on the 
right hand side is restricted to p. If there is a non-trivial homomorphism A : F — > G, it will 
induce a corresponding homomorphism A : F_ — > G on the Lie algebras and we can choose 
the Wm(y°) to represent it. 

Next, at yo, consider the (15. 5115. 6p for m. Eliminating L^W M , and noting that alone 
depends on x, we must have (a) \W m W n ] = f p mn W p and (b) [<3> n , Wy = 0. Note that this 
implies that the residual gauge group is reduced to those elements of G which commute with 
W m i.e. to the centralizer of X(F) C G. The gauge transformations are already restricted to 
be functions of x alone. 

Considering the Jacobi identity for $j, it follows that = d m kn Q m must hold 

for some d's. This has the same form as the condition (a) on the W's. Hence d m kl = /™ z is 
obviously a solution. 

In fact, it is a result that S— invariant connections, when they exist, are in on-to-one cor- 
respondence with homomorphisms of the groups A : S — > G and can be expressed as 
A(x,y) = Aidx 1 + §(x) n i*(QMc) n (y) where A is a connection on B with the gauge group 
reduced to the centralizer of X(F) in G (i.e. group of all elements of G which commute with 
the image of F in G under the homomorphism A). Furthermore the Higgs fields have to 
satisfy the constraints: $n] = f P mn^p^ 

What about other invariant fields, such as vector fields in the adjoint of the gauge group (eg 
the triad fields)? Now the invariance condition (15.51) will change and also the corresponding 
gauge transformations of the field. For E^d^, we will have L^ m E = [E, W m ],E 9 = g~ x Eg. 
Now use the projections X m of the left invariant vector fields X m on S (these generate the 
right action and are dual to the 0mc) an d write: E" := ^"X™, L^ m X n = 0. The invariance 
condition then becomes {u k a X*){L u ^ n - W m )) = and {uj k a X^){{^ n f - g~ l ^ n g) = 0. 
The gauge transformations imply \E fn transforms by the adjoint representation and for m = 
m, the invariance condition implies: [$ n , W^j = 0. Exactly as before, the \l/ n must satisfy 
constraints analogous to the $ n 's. Similar logic will hold for other tensor fields. 

As a very simple illustration, consider the case of static magnetic field in three dimensions 
invariant under translations along the z-axix. We want to obtain the form of the vector 
potential A*. In this case, the gauge group G = U(l) and the symmetry group S = E 
which acts on IR 3 by translations. This action is free and therefore F = {e}. We have 
]R 3 ~ ]R 2 x K. The Maurer-Cartan form on IR is just dz. The map from R — > U(l) is given by 
a single 'Higgs' scalar, $(x,y). The symmetric connection is then given by Ai(x, y, z)dx l = 

4 The classification of connections invariant under some group of automorphisms of appropriate bundles is 
given by generalised Wang theorem. There are many mathematical fine prints in the above discussion which 
should be seen in the references in [23]. 
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A x (x,y)dx + A y (x,y)dy + <p(x,y)dz. This just says that all the three components of the 
vector potential depend only on (x,y). Note that this is a statement in the gauge where 
'W has been set to zero. This implies that the magnetic field is also independent of z. 
Note that since A z := <p(x,y) ^ 0, the magnetic field could be along any fixed direction. 

Exercise: Work out spherically symmetric Yang-Mills fields in three dimensions. Now G = 
SU(2), S = SO(3), F = 17(1), M 3 ~ M+ x S 2 . 

Further examples may be seen in [23J. 



5.2 Schrodinger and Polymer Quantization 

We illustrate inequivalent quantization as well as the GNS procedure in a simple example. 



5.2.1 The Weyl-Heisenberg C*- Algebra 



Consider the usual Schrodinger quantization of a single degree of freedom. We have the 
usual Hilbert space L 2 (M, dx), on which are defined the self-adjoint operators x,p, satisfying 
the canonical commutation relations: [x,p] = ih. 



Define the corresponding unitary operators: 

U(a) := e iax , V(p) := e ih ~^ p , a,(3e 



(5.9) 



a) 



U(-a) = Uia)- 1 , V\(3) = e'^P = = V(P)~ X 



Using the BCH formula, 



it follows, 



W{z) 



It follows, 



e A . e B = e A+B+l[A,B] + ^[A,[A,B]]-j2[B,[A,B]]- 



U(a)U(a') = U(a + a') , V(P)V(P') = V(/3 + /?') , U(a)V(/3) = e - ia ^V(/3)U)(a) 

(5.10) 

Define, for z := (a + i/3)/y2 G C, the unitary operator, 



e l -U(a)V((3) 
exp 



z + z z — z p 
% < — -^x — % 



y/2 



V2 h 
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x — ip/h _x + ip/K^y 
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V2 J 
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1 / , .P 
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e -4- e i2at e i2a where, 
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(5.11) 
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e^ zlZ2 'W(z 1 + z 2 ) 
W(-z) = W(z)- 1 



and, 



(5.13) 
(5.14) 
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Taking finite linear combinations of products of the unitary operators W(z), we get an 
algebra called the Weyl-Heisenberg algebra, W. This is *-algebra due to the Hermitian dagger 
defined for the operators. The unitary operators W(z) are bounded and so are polynomials in 
them. With respect to the operator norm (which satisfies 1 1 \ \ = \\A\\,\\A^A\\ = ||A|| 2 ),The 
Weyl-Heisenberg algebra is a C*-algebra. Notice that the W C*-algebra is non- commutative 
and has two commutative sub-algebras, namely those generated by the elements, W(^Q) 

and W(^j£-) respectively. 

Thus at this stage we have constructed a C*-algebra of bounded operators on the specific 
Hilbert space. We will now define a positive linear functional on the C* algebra, W, con- 
struct a unitary representation of the algebra and show its equivalence to that provided by 
the Schrodinger quantization. The same procedure will then be used to construct another 
representation, the Polymer Representation, of the same algebra. 



5.2.2 Re-construction of the Schrodinger Representation 

In the Hilbert space, consider the wavefunction, (x\0) := ipo(x) := that 
(0 1 0) := j dx\ip (x)\ 2 = I. Following results hold: 

a\0) = -j= (x - Vo(ar) = (5.15) 

[W(z)4> ](x) = e^[U(a)V((3)^ ](x) = e 5 [V(^ ] (x) 

= e^e iax ^ (x + (3) (5.16) 

= (x\W(z)\0) = e-^(x\e lzat \0) := (x\z) (5.17) 
.-. J dx%{x) \W(z)if> ] (x) = ir- 1/2 J dx e'Tj 



,2 . Q+/3) 2 
tax— k 



M 2 

e~— ■= (0\W(z)\0) = (0\z) (5.18) 



Define a linear functional fisch, on the Weyl-Heisenberg algebra 

fisch fec^( Zl ) := ^cAschW*)) : = J>(0|W^)|0> := ^Qe^l 2 / 2 . 

\ i / i i i 

The first definition ensures linearity and the third one completes the definition. The sec- 
ond definition (notational) makes it obvious that fis c h is a positive linear functional since 
Q S ch(A*A) = (0|AU|0) = ||A|0)|| 2 > 0,V A : = CiW(zi) e the C*-algebra. The equality 
holds only if A\0) = Ci\zi) = 0. Since \zi) states are linearly independent, there are no 
non-trivial states in the Hilbert space which satisfy A\0) = 0. The positive linear functional 
then defines an inner product on the algebra by, 

(W(z),W(z')) := n Sch {W{ Z yW(z')) = (0\W{-z)W(z')\0) = (0\z' - z) . 

and extended by linearity to the algebra. This turns the algebra into a Hilbert space (distinct 
from the original Hilbert space). 
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Next, define operators, W(z) acting on the algebra, by, 

Im(zz') 

W{z)[W(z)} := W(z)W(z') = e^^W{z + z') 

and extended to the algebra by linearity. Similarly, one defines an operator for each element 
of the algebra, in an obvious manner. 

Exercise: Show that {W(z)Y = W(—z). 
This implies that W(z) are unitary operators. 

That W(z) —> W(z) provides a homomorphism of the W algebra is obvious from the action 
of the operators. Thus the algebra with the inner product defined, carries a representation 
of itself in which the W(z) are represented unitarily. 

Consider general matrix elements of the operators W(z): 



W(z 1 ), W(z) \W{z 2 )\j = (0\W(-z 1 )W(z)W(z 2 )\0) 

= e 2 e 2 (5.19) 

Observe that for z = a or z = i/3, the above matrix elements are continuous in a, (3 respec- 
tively. General matrix elements are obtained from finite combinations of these and hence the 
W(a) and W{if3) are both weakly continuous families of unitary operators. Actually, these 
are also strongly continuous families i.e. w.r.t. vector space norm. The strong continuity 
can be checked by evaluating the norm 1 1 ("W^(^) — W / (0))[W r (z')]|| and checking the limits for 
z = a/\/2 and z = i/3/ \/2. 

This allows us to define two self-adjoint operators (on the algebra) as, 

X := to^W-' , P := H (5.20) 

«->o ia /3^o i/3 

X,P = ihl (5.21) 

The commutator can be evaluated directly using the definitions for a, (3 ^ and using the 
existence of the limits guaranteed by strong continuity. 



5.2.3 Another positive linear functional and the Polymer repre- 
sentation: 

Now view the Weyl-Heisenberg algebra defined above as an abstract structure i.e. an algebra 
generated by elements W(z),z G C, obeying the relations (I5.13|5.14p with a norm defined 
by ||iy(z)|| = 1 V z G C and extended by linearity. Define a linear functional by, 
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This is positive because 

Y,c*c 3 e-^ z ^n moy {w{z 3 - Zi )) 

E W + E CtCj^-** 6 ™ ( 5 - 23 ) 

In the sum, only those pairs (2, j) which have the same /3, contribute. Group together all 
the terms whose are equal (eg {zi, Z2, ■ ■ ■ , z m }, {z m+ i, . . . z m+n } . . .), and consider one such 
group at a time. In each such group, the phases in the second term can be absorbed in the 
Cj's (since the (3 is common) and then combining with the first term gives | '^2 i {Cie' ll3a ^ 2 }\ 2 
which is non-negative. This completes the proof. Note that positive linear functional must 
evaluate to 1 on the identity element of the algebra (namely, I := VK(0)) and therefore we 
cannot interchange •<->■ a in the defining condition in (I5.22p . 

The fipoiy defines a degenerate inner product on the algebra, 

(W{z'),W{z)) := n Poly (W(z'yW(z)) = e- 1J ^n PolY (W(z-z')) 

and extended by linearity. The elements whose norm w.r.t. this degenerate inner product 
is zero, forms a closed subspace Af, of the algebra and consists of elements of the form 
X ■■= Y.i C i W (.^f)^ 6 Msuch that J2i C i e ~^ ai/2 = °- Elements of Af also satisfy the 
property: Q Po \ y (Ax) = V A e W, which is useful in the exercise below. The quotient 
space, W/Af, is an inner product space and its Cauchy completion defines a Hilbert space 
of the Polymer Representation. 

Exercise: Let A denote a general element of the algebra and x an element of Af '. Define 
[A] := {B E W/B = A + X }- Define ([A], [B]} := (A, B) and := [W{A]}. 

Show that these definition are well defined and conclude that W/Af provides a unitary 
representation of the quotient algebra. From now on, we refer to the quotient representation 
without being explicit about it. 

Observe that flp \ y is continuous in a and discontinuous in /3. This directly implies that 
W (ij3 / cannot be weakly continuous and therefore we cannot define the analogue of P. 
This follows by noting that (W(z'), W(i/3/V2){W(z')}) ~ tt Poly (W(i/3/V2 + z' - z')). The 
weak continuity (actually also strong continuity) in a however allows the definition of X 
self-adjoint operator. It remains to make the representation explicit. 

In both the cases above, with the Schrodinger and the polymer functionals, we constructed a 
representation of the W in which the W (z) are represented by unitary operators. This is the 
Gelfand-Naimark-Segal (GNS) construction. In the Schrodinger case, we obtained the X, P 
operators satisfying the canonical commutation relations. In the Polymer case we obtained 
only X. In both cases we get the following relations, directly by applying the definitions: 

W ] X W (ip/Vt) = X-pi (5.24) 

XW (ifi/V2 \ - W (i/3/V2) X = -0W {i/3/V^j (5.25) 

From these relations follow an important result. Note that X is a self adjoint operator and 
therefore its spectrum is real. What can one say about its eigenvectors? The above relation 



^pioy E owe*)}* E i C i W te)} 
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implies that if \x) is an eigenvector of X with eigenvalues x, then W(ip/y/2)\x) is also an 
eigenvector with eigenvalues (x + j3). Hence, either every real number is an eigenvalue or 
none is. 

Suppose every x is an eigenvalue. Then we have the orthogonality relation (x\x') = 8 XjX > 
- the Kronecker 8. Taking expectation value of the second equation above, it follows that 
f3f(x,/3) = 0,Vx,/3 E R, where f(x,/3) := (x\W(ip/\/2)\x). This implies that f(x,/3) is zero 
for (3 7^ and f(x, 0) = 1 directly from the definition. Thus f(x, (3) cannot be continuous at 
— 0, for any x. This means W{i0 / \/2) cannot be weakly continuous at = 0. 

This also means that if W(i/3/y/2) is weakly continuous (as for the Schrodinger represen- 
tation), then X cannot have any eigenvector. Each x G R is a generalised eigenvalue and 
hence, in the formal notation, (x\x') = 5(x — x') - the Dirac 5-function. 

Thus, in the Schrodinger representation, X necessarily has only generalised eigenvalues, while 
in the polymer representation, it could have proper eigenvalues, but generalized eigenvalues 
is not ruled out. 

However, in the polymer representation, we note: 

r W{ip/V2), W(i0'/V2)) = fipoi y (w(-ip/V2)W(ip' 



W(i0/V2), x{w(i/3'/V2)} 
(X- XS)\W(ip/y/2) 



P' 





1 " 


= lim 




a->0 


ia . 



W{ 







a+ip > 

for A 



= 5 Pt p ; (5.26) 
(5.27) 

1 + ia\)] W(if3/V2) 
P . (5.2J 



which show explicitly that W(—ift/y2) is a normalized eigenvector of X with eigenvalue (3, 
for every G R. This means that the Polymer Hilbert space is non-separable. 

This concludes the illustration of the GNS construction of representations of C* algebras. 
In the next sub-section we see the analogue of the spin network construction. 



5.2.4 Polymer representation via 'spin networks' 

We begin by introducing 'graphs' in a '0-dimensional manifold', define 'holonomies' and 'spin 
network functions', define an inner product and densely defined operators. More details may 
be seen in [21]. 

1. Graphs, holonomies, cylindrical functions: Any countable set of real numbers, {xi} 
represents a graph and is denoted by Note that the 'points' Xj on the real line 
correspond to edges. Associated to each edge, Xj, we define a point holonomy, e lkXj , 
fceM plays the role of a connection. For each graph 7, define complex valued functions, 
/ 7 (A;) := Y^,j ffi %kx i ■ Let Cyl 7 denote the vector space of f^ik). Elements of this vector 
space are said to be functions cylindrical with respect to the graph 7. Let Cyl := 
CyLy, where the sum is over all possible graphs. 

Thus a general element of Cyl is a function of k expressible as a countable linear 
combination of the elementary functions f Xj (k) := e J 's. 

5 For precise technical conditions, please refer to [24] 
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2. Inner Product on Cyl: Define 



(f,g)poi y ■= lim — / dkf*(k)g(k) 

L— >oa ZLi J _l 

which for elementary functions f x (k),f Xl (k) gives 5 XtXr Introducing the notation, 
\xj) -H- e lkXj this is expressed as (xj,x{) := S X tXl , V Xj, Xi G R. Cauchy completion of 
Cyl with respect to this inner product defines a Hilbert space, "Hp i y . 

3. Action ofW: Define operators W(z) on Cyl by, 

W(^)f{k) := e-^ 2 e^ k f(k-a) , V / G Cyl, V a, f3 G R 

These are densely defined and can be extended to bounded unitary operators on "Hp i y . 
It is easily verified that this provides a representation of the abstract algebra defined 
in (15U31) . 

Exercise: Show that the 1-parameter families of unitary operators, W(a/y/2) and 
W(i/3/y/2) are weakly continuous at a = and weakly discontinuous at f3 — respec- 
tively. This implies that while X can be defined from the W(a/y2) family, there is 
no corresponding operator from the second family. Thus the holonomies - h x - (k) are 
well defined but not the connection - k itself. 



5.2.5 Harmonic Oscillator in the polymer representation 

So we see two distinct representations of the same abstract algebra with the polymer rep- 
resentation being similar to the LQG representation. Are there observable quantities which 
would reveal which representation occurs in nature? 

An obvious candidate is to consider the dynamics of the Harmonic oscillator, with the clas- 
sical Hamiltonian, H(x,p) := p 2 /(2m) + mu 2 x 2 /2, {x,p} = 1. In the quantum theory, the 
x,p are expected to be replaced by the corresponding operators. However, in the polymer 
representation, there is no p\ So in proposing the quantum Hamiltonian we need to intro- 
duce a scale, no and define p 2 := {2 — W(ifio/y/2) + W(—in /\/2)}/fil. We could of course 
define p first and then take its square. This is a quantization ambiguity not too important 
for our purposes here [23] . In the Schrodinger representation, we could use exactly the same 
definition, work out quantities of interest eg spectrum and then take the limit (1q — > 0. In 
the polymer representation, we cannot take this limit. 

Consider the eigenvalue equation for the Hamiltonian: H\ip) = E\ip). Writing \if>) := 

e countable set \ x ) > an d noting that the unitary operators in p 2 shift the \x) — >■ \x±/i ), 
the eigenvalue equation becomes a difference equation, involving ip(x), ip(x±fi ). This means 
that ip(x) with x in a lattice L xo := {x = x + HqN, ]VgZ} constitute a solution while those 
belonging to different lattices are unrelated. Span of the vectors in any lattice form a sep- 
arable subspace of the polymer Hilbert space. The spectrum can be determined for each 
lattice independentally. This is analyzed in detail in [23] • Suffice it to say that the spectrum 
differs from that in the Schrodinger representation by terms down by powers of /xo/d. Here, 
d:= y/KJ mw is the length scale defined by the system while fiQ is the length scale introduced 
by the approximation for the momentum operator. For the physical systems modelled well 
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by an oscillator (eg for vibrational spectra of molecules), [io/d is extremely small and so 
Schrodinger vs polymer representation cannot be resolved by observations. 

Additional comments may be seen in [24] . 



5.3 Inverse Triad Operator(s) 



As noted before, the discrete nature of the spectrum of the triad operator implies that its 
inverse is not densely defined. Consequently the counterpart of the classical function p^ 1 
needs to be defined indirectly, by a suitable prescription. Being a prescription, it introduces 
quantization ambiguities. We will consider a prescription which is sufficiently general. 

We aim to define an operator sgn(p)|p| _1 . Introduce the following notation: n l is a unit, 
3-dimensional vector and Tj are anti-hermitian generators of 577(2) in the J th representation, 
satisfying 

[t^Tj] = e ijk r k , Tr J (r l r j ) = ~J{J + 1)(2J + 1)<% := - jVjtfy 

For f a (p) : = /i a |p|- Q / 2 , and define g a (p) := J p f~ 1 (x)dx = n~ l ?>ga{p)\p\ 1+a / 2 (l + a/2)" 1 . 
For a = we have fo = fio while for a = 1 we have f\ = ^i\p\~ 1 ^ 2 and we choose /ii : = 



£py7v^3/4. These recover the /i and the /z schemes. We will suppress the a label. Define 
hf := e AraV i/( 0! .p) c . This is matrix of order (2J+ 1). Classically, the following is true. 

¥\f^j ■ d I CJ I 

hf{h~^ 1 , |g|'Go)} = — — (AnVj) Z|(7|' _1 sgn(p) where we used /— — = sgn(p). (5.29) 



| p |(i-l)(l+«/2) = sgn(p) 



sgn(p) 



(^(l + a/^y-^J 1 



x [Trj((nW t )h f {h]\\g\ 1 })} (5.30) 



^(/i a (l + a/2))A0l 



x [^((^^{^Mp]^/ 2 )})] (5.31) 

Thus we have a classical expression for |p|(' _1 )( 1+a / 2 ) which has four ambiguity parameters: 
a, J, I, A. J is a positive half integer, < I < l,a > —2. The special cases would be: 
(a) /i — scheme: a = 0, A = 1; (b) improved scheme: a = 1, A = 1,7 = 1/2 and some 
special values of / explored; (c) Madhavan scheme: similar to the improved scheme except 
A is correlated with the fiducial volume Vq (more on this later). We could choose I — 1 = 
(1 + a/2)~ 1 to define inverse triad, but we will postpone such choices. 

The corresponding quantum operator is obtained by replacing the Poisson bracket by —iftr 1 
times the commutator. The —i is combined with n ■ r to make the generators Hermitian and 

the hr x combines with k to replace k by £p. The commutator is expanded as: hflhj 1 , \g\ ] = 
I-bl - hf\g\ hf 1 . 

Observe that n ■ (—it) can be diagonalised with diagonal elements being — J, — J + 1, ... J — 
1, J. So the hf becomes the diagonal matrix e^ lX ^ aC ^ J ' J ~ 1 ' "• - J+1 - J ), So the commutator 
terms are diagonal matrices. 
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The computations simplify if we label the basis states by 

So that the /i/ shifts the v a labels simply as, 

e^ c \v a ) = \v a + 2k\) y e i/c/ » = \v+ 1) . 

With these, acting on a basis state \v a ), the Trj evaluates to, 

i J 



[TrA---}] = j^k{\v a \ l -\v a -2kX\ 1 } 
^ ' k=-j 

= ( ~^r~ ) 5^ ^l^a + 2fcA|' and defining ?;„ := 2JXq a , 

= (l^y#-£ k \Jq a + k\ l := ^^y 8gn ( fc ) GjiJ ( fc ) (5.33) 
J 

Gj,,(g a ) := sgn(g a ) ^ k\Jq a + k\ l (5.34) 

fe=-j 

The eigenvalues of \p\( l ~ 1 )( 1 + a / 2 ) are then given by 

\ p \<MW+a/2)\ Va ) := S gn(v a )Aj, l>a (v a )\v a ) (5.35) 

AjM = (^y(i±jf^^ [GjM 

The first bracket takes care of the dimensions and the remaining factors are dimensionless. 
It remains to calculate the last square bracket which is a A, a independent, universal function 
of its argument and depends only on J, I. 

From its definition, it is easy to see that Gjj(0) = and Gj t i(—q a ) = Gj t i(q a ), the sgn(g a ) 
factor is crucial for this. Thus it suffices to consider only q a > 0. 

j 

G jt i(q a > o) = k \k + Jq a \ l 

k=-J 

J 

= Yl i( k+ J ^)\k + Jq Q \ l ~ Jq a \k + Jq a \ 1 } 

k=-J 
J(q a +1) 

= Yl {sm(k)\k\ l+1 - Jq a \k\ 1 } (5-36) 

k=J(q a — 1) 
J(q a + 1) 

Gj,i{q a >l) = {\ k \ l+l ~ J V<M 1 } (5-37) 

k=J(q a — 1) 

0_ k=J(q a + l) 

G jil (0<q a <l) = {\k\ l+1 + Jq a \k\ 1 } + E {\k\ l+1 -Jqa\k\ l }(5.38) 

k=J(q a -l) 0+ 



6 The functions f a (p) are taken to be dimensionless. This makes the ii a to have dimensions of £p and 
9a{p) to have dimensions of £p. The v a is defined to be dimensionless. 
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In the second step, we have shifted k — > k — Jq a . k is no longer integral but still changes 
in steps of 1. Clearly, for q a > 1, k is positive and the sgn as well as the absolute value are 
redundant (the k = term for q a = 1 gives zero and hence the sum is confined to positive k 
only), as in (I5.37p . For q a < 1, the sum splits in two groups as in (I5.38p . and the 0± denote 
the respective limits on the values of k which must match with the other limits and shift in 
steps of 1. 

We will not simplify/approximate this further but consider the special cases (i) a = 0, A = 1 
and (ii) a = 1, J = 1/2,1 = 2/3. 



(i) a = 0, A = 1: (eigenvalues of \p\ l x ) 

fo(p) = t*o , 9o(p) = Vo ls g n (p)\p\ , ^o = ^o ls 6 n (^)l^l 

Aj A o(%) = 0^) (/A0)- 1 G J XI^|/(2J/i )) 



A i i M) 



2'2 



-1/2 



ii 



+ 1 



1/2 



ii 



l/2^ 



(5.39) 
(5.40) 
(5.41) 



(ii) a = 1, J = 1/2: (eigenvalues of \p\i {l 1} ) 



hip) = 1/2 



f'r.- M 1 ^ 3 -^ , 9n{p) = ^g<P)\P? / \m- 1 



6 2 



^y / Vr 1 (3/2)- 1 sgn( yU )|/i| 3 / 2 = ^sgn(/i)H 3 / 2 



Ai (1 K 



^1743 
3 2 



i-i 



A (lATi)- 1 G? 1/sy (h|/(2iA)) 



The G can be computed directly from the (15.341) . for v\ > 0, as 



l+l 



ii 
A 



+ 1 



For 1 = 2/3 the operator becomes \p\ 1//2 and the eigenvalue becomes, 



Ai 2Avt) 

2'3' 



6 



2 \ -1/2 



3 fK\ 3 



4 \ A 



T + 1 



(5.42) 

(5.43) 
(5.44) 

(5.45) 



(5.46) 



For large v± and A = 1, this matches with the eigenvalue given by [20)1 . The difference arises 
because the APS prescription takes the a = expression and replaces /i by p, in equation 



For large J the sum can be approximated using, 

f 1 1 N { ' V 1 N 

and applying it to the sums in the definition of the Gji. 



N 



r+l 



r + 1 



7 Actually, for large volume, the leading term is independent of A. The sub-leading (correction) terms, do 
depend on A. 
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5.4 Inhomogeneous Lattice Models 



In the main body, we focused on quantization of symmetry reduced models which are based 
on a homogeneous and isotropic background. From a perspective of the full theory, this 
background presumably corresponds to a state of full theory. Generic states of the full 
theory would be inhomogeneous. One way in which symmetric states of the full theory have 
been understood in the LQC context is that the symmetric states are those distributions 
in Cyl* of the full theory which have support on the invariant connections [23]. In the 
same spirit, we may stipulate certain kinds of inhomogeneous states as those distributions 
which have support on certain form of 'inhomogeneous connections'. Specific models can 
be then constructed using similar strategies as used in LQC constructions. Such models 
can shed some light on how homogeneous and isotropic models could be viewed from an 
inhomogeneous perspective. These so-called lattices models are briefly summarised below. 
Details should be seen in [16]. 

For definiteness, let us continue to work with homogeneous (not necessarily isotropic), spa- 
tially flat, diagonalised model with a fiducial cell of co-moving volume Vq as before. The 
spatial isometries provide directions (of the Killing vectors) and the fiducial metric provides 
coordinates as background structures which are to be kept fixed. Using these background 
structures, Bojowald constructs another model as follows. 

Choose a cubical lattice (say) aligned with the isometry directions and with a spacing £ := 
(Vo/iV) 1 / 3 . Let the vertices of the lattices be denoted by v and the three oriented links be 
denoted by ej^(t) := v + tej, t e [0, £ ) and ej is the unit vector in the I th direction. For 
future reference, let Si $ denote the elementary surface perpendicular to the elementary link 
ejj and passing through its mid-point. 

Restrict the connections and triad variables to be of the form, 

4(x) := ~k I (x)Sl I / a , E?{x) := p\x)5f . (5.47) 

These are the local versions of the diagonalised homogeneous models. The diagonal form of 
the connections implies that the holonomies - path ordered exponentials - become ordinary 
exponentials of line integrals. 

The kjjp 1 are further taken to be spatially periodic with period V^ 3 , 
kj(x) = ^2ki{m)e m - 3 , p\x) = 5^(m)e ifll,5f , m = 2nV~ 1/3 n, n 6 Z 3 . (5.48) 

in rh 

The Poisson brackets between the connection and the triad lead to, 

{k I (x),p J (y)} = K^5j5 3 (x,y) . =4> {kj(m),p J (m')} = K,<yV ~ 1 8j8 3 (m, -ml) . (5.49) 

In loop quantization, basic variables of the model will be holonomies of the lattice connection 
along the three elementary links at each vertex and the three fluxes of the lattice triad 
variables along the elementary surfaces through the mid-points of the elementary links and 
perpendicular to the link. As noted above, these holonomies will be ordinary exponentials 
thanks to the diagonal form of the connection. 
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The line integrals of the connection along elementary links of the lattice are given by, 
Zj,* := / dt h(ei(t)) = f dtJ2~ k i(™) el ™'*'' v<t) = himy™-* j ° dte**** 1 
= J2h(m) [2 e «fl-V"/W2 (" Sin( ^ o/2) ") j where, mj := m • e, . (5.50) 

~ kj(v)£Q •/ m/£ ^ 1 dominates the sum. 
/i/^ := e2 X/ .* « e 5*/Wo ;= e 5*/W (elementary holonomies) (5.51) 

Likewise, the fluxes of the lattice triad along elementary surfaces are given by, 

Jl ■= \ ^p J (m)e^ = Yv J {jn)e^e im ^l 2 \ e ltmK dt / e ltmL dt 

JSjj A ^ J-e /2 J-io/2 

= Ue^e^ 2 ^in(m^V2)sm(m £ f /2)\| ^ 

rh 

pi(v)£% (elementary fluxes) (5.53) 

The J, K, L indices are chosen such that €jkl = 1- This takes care of the orientations. There 
are no smearing functions above because the p 1 variables are (U(l)) gauge invariant thanks 
to diagonalised form. 

These variables satisfy the Poisson brackets, 

j rh .(^> ) sm(m I l /2)sm(m K l /2)sm(m L £ /2) 



= K'ySj \xto(v ~ ^)] = Kl&ih* ( 5 - 54 ) 



m 



The square brackets above is the characteristic function of width £q and centered at (v — v') 
which is just the Kronecker delta. 

The kinematical Hilbert space is then described in terms of the flux representation as: 

H\ •••» ^i,v , •••) = ( -^-m,v\ I in# , •••) , in$ e z (5.55) 

^/,^| • • • , m,v , • • •) = I • • • , Hi,* + 1 , • • •) (5.56) 

The flux eigenvalues are in integer steps of j£p/2 because the elementary holonomies suffice 
to separate the lattice connections (periodic) and thus only their integer powers appear. 

Subsequent steps are similar to what is done in the homogeneous models. In particular, the 
volume corresponding to the cell with N 3 lattice sites, can be expressed as 

V V 

- E (5-57) 

V 

leading to the corresponding operator expression. 
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There is no diffeomorphism constraint since the background coordinates are fixed in defin- 
ing the lattice, the SU(2) gauge invariance is first reduced to the U(l) 3 due to restriction 
to diagonal connection and triad and by the form of these variables, the k^p 1 are gauge 
invariant variables. Hamiltonian constraint remains as in the case of homogeneous models. 
In essence, we have N 3 'homogeneous models' (labelled by the Fourier label m) at the level 
of basic variables and the kinematical Hilbert space. The role that Vo played in the homo- 
geneous model is now played by £ . The inhomogeneity is reflected by basis states having 
different values of variables. 

How do we relate this set-up to the isotropic one discussed before? 
Observe that a generic basis state in the lattice model will be, 

VWife] = niftif)" 1 '* := ( k A x )\ ■ ■ ■ > , ■■■) , k A x ) := kj(x)£ . (5.58) 

I,v 

~ -I /O 

If we choose k(x) := c := V c V x,I, then the basis function becomes a function of a 
single variable c (which is independent of x), and is of the form: 

M<>) = > fi'=N- 1/3 ^2fi Ifi GQ; (5.59) 

r,v 

which can be viewed as a basis element of Cyli SOtr0 pic- Thus we can define a map 7r : 

Cyllattice ~ ^ Cylisotropio 

7T : I , ...) -> \fi) & (c\fi) := (kj(x)\ , . . .) \~ k(x)= - c 

with /I := A^ 1 / 3 ^/^ (5.60) 

I,v 

Note that the image of 7r-map is a separable subspace of Cyl iso t r o P ic, spanned by e Q. 

Clearly we cannot uniquely identify a cylindrical state of the lattice model, given a cylindrical 
state of the isotropic model. However, we can define a map o : Cytotropic — > Cyl* attice , 
a : |//) — > (jl G R, such that, 

(//|..., •••) = (/x|7r(| i//,* , ...)) = , i/ := iV" 1/3 "i,v ■ (5.61) 

In the second equality, we have used the inner product of the isotropic model. This map 
embeds cylindrical states of the isotropic model into the distributional states of the lattice 
mode^. 

Now, we have Operators A* acting on Cyl{ attice corresponding to operators A acting on the 
Cyli a ttice, defined in the usual manner. Those of these operators which act invariantly on the 
image of a in Cyl* attice , can be identified with operators of the isotropic model. For these 
operators, we can define isotropic via the equation: a{A motvop i C \ii)) := (cr|/i))A* attice . Since 
we have embedded isotropic states in the distributions of the lattice model and also have 
correspondence between operators, matrix elements computed in the isotropic model can be 
understood as actions of lattice distributions on lattice cylindrical states. 

8 Notice that /x, v defined above are rationals with a common denominator TV 1 / 3 . Therefore the distribu- 
tions are non-trivial only for /igQ with the same denominator. 
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Consider an operator Ai att ice on Cyli at tice- This defines an operator A* atticc on Cyl* attice : 
(^lattice^l • • • » u i,v , •••) : = (01 {^lattice I "if, •••)}• If for every distribution (0| = 
(fx\ =: a (\ (/,)), the operator -A* attice gives another distribution =: cr\ [/,'), then we get an 
operator on Cytotropic : 

^■isotropic | A*) • IA*)- 

It is easy to see that the multiplicative operators on Cyli at tice, give multiplicative operators 
on Cytotropic- For example, taking Attice = hj^i, the lattice state | . . . , , . . .) will have 
the vjfi incremented by 1. The action of the {p\ will give S^ v+ i. This can be understood as 
the action of (/x — 1| on the original lattice state. Thus (A* attice /x| = (^ ~~ 1| which implies 
the a multiplicative action A isotropic |/x) := |/x — 1). 

For elementary flux operators, little more work is needed. For example, action of T\ on a 
basis state, | . . . , vj$i , . . .) is zero unless vj$i ^ for some .1 = 1 and at some tf = v i.e. 

^Wj,v>) = \lfipVj,v>5j8^v\vj$,) or 



$j$v,v'$h,v ' v := N x l z vjtf However, (5.62) 



2 

(fjf\vi,v>) = (fjf\u Iii?l ) V (/x'|, v, v" . (5.63) 

Thus {Pi)* cannot act invariantly on the image of a in Cyl* attice . It is clear though that if 
we sum the elementary flux operators (with the directional index J) over all the lattice sites 
(this is a finite sum due to the cell), then the sum will act invariantly. By averaging over 
the directions as well, we can construct an operator corresponding to a 'flux' operator on 
Cytotropic- In equations, 



^lattice 



N' 1/3 J2^ > ^tice := ^lattice =► (5-64) 

(A*| {Platticel • • • , Vjf , . . .)} = -74A r_1/3 ^2 • • • ' V W ' • • •) 

J,v 

I,v 

The last expression is exactly the matrix element of the p operator defined in the isotropic 
model. As noted in the footnote, the identification of the matrix elements is restricted to 

This completes our summary of the lattice model and how isotropic model is 'embedded' in 
the lattice model. This 'embedding' refers to embedding of the particular separable subspace 

Of Cytotropic- 
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